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^ ■ Abstract 

Q-i! The simultaneous investigation of the pion electromagnetic form factor in the space- and time- 

p ,' 

Q I like regions within a light-front model allows one to address the issue of non-valence components 



of the pion and photon wave functions. Our relativistic approach is based on a microscopic vector 



^ . meson dominance (VMD) model for the dressed vertex where a photon decays in a quark-antiquark 

\ pair, and on a simple parametrization for the emission or absorption of a pion by a quark. The 

results show an excellent agreement in the space like region up to -10 {GeV/c)'^, while in time-like 
region the model produces reasonable results up to 10 iGeV/c)'^. 
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I. INTRODUCTION 



e.g. 
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Electroweak properties are widely used as an important source of information on the 
structure of hadrons. In particular within the framework of light-front (LF) djTiamics 

, Sl^a large number of papers has been devoted to the study of nuclei and hadrons (see 

\iu innn 

19|, [201, [211, 1221 , just to give a partial account 



^ large numoer oipapei 
QBIallO, U, 12, 13, [l4 
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of previous works with a finite number of constituents). 

The LF dynamics allows one to exploit the intuitive language of the Fock space. Indeed 
the Fock-space language is particularly meaningful within LF dynamics, since: "The sim- 
plicity of the light-cone Fock representation as compared to that in equal-time quantization 
is directly linked to the fact that the physical vacuum state has a much simpler structure on 
the light-cone because the Fock vacuum is an exact eigenstate of the full Hamiltonian." P| 

Another basic motivation for choosing the LF dynamics is represented by the striking 
feature that the Fock decomposition is stable under LF boosts, since they are of kinematical 
nature and therefore do not change the number of particles, i.e., are diagonal in the Fock 
space. 

Therefore the LF dynamics is a suitable framework for the investigation of the Fock 
expansion for mesons and baryons, viz 



\meson) = \qq) + \qqqq) + \qq g) + 

\haryon) = \qqq) + \qqq qq) + \qqq g) + 



(1) 



In particular, within the LF dynamics the e 



ectromagnetic form factor of the pion has 
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2l|). Indeed the 



been the object of many papers (see, e.g., Refs. 3, 9, 
pion electromagnetic form factor yields a simple tool for the investigation of pion and photon 
microscopic structure in terms of hadronic constituents. In what follows we will present an 
approach to investigate in a common framework the pion and photon vertex functions, with 
the perspective of an extension of our approach to the nucleon. The intuitive language of 
the Fock space will be widely exploited to analyze the above mentioned vertex functions. 

Aim of this work is to give a unified description of the electromagnetic form factor of the 
pion, both in the space-like (SL) and in the time-like (TL) regions, taking into account the 
complexities related to the pion and photon vertex functions, both in the valence and in the 
nonvalence sectors, as well as the fermionic nature of the constituents. A first presentation 
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of our approach was given in Ref. 23 1. 

The choice of the reference frame where the form factor analyses are carried out has a 
fundamental role, as shown in previous works in the space-like region 

Q Q, Q, Q 

and 

in the time-like one [itI]. For a unified description of TL and SL form factors, a reference 
frame is needed where the plus component of the momentum transfer, g+ = g° + g^, is 
different from zero (otherwise, = q'^q' — q]_ cannot be positive). As a matter of fact, a 
reference frame where g"^ 7^ allows one to a naly ze, in a common framework 2^, the pair 



production process (Z-diagram contribution) 18|, i.e. the effect of multiquark propagation, 



as well as the ultrarelativistic effect of the so-called instantaneous contributions and the 



hadronic components of the photon wave function 

In Ref. j26[ it was shown that, within the Hamiltonian LF dynamics (HLFD), a Poincare 
covariant and conserved current operator can be obtained from the matrix elements of the 
free current, evaluated in the Breit reference frame, where the initial and the final total 
_a of ...en, d.ec.ed alon. .Ke .p,„ ,uant.a.o„ . Following Ref. Q. 
we calculate the pion form factor in a reference frame where q 1 = and 0+ > 0. 

n 

Our starting point is the Mandelstam formula [23] . To construct a bridge toward the 
Hamiltonian language, the hadron vertex functions will be connected to the LF wave function 
of the valence component of the hadron state. Furthermore, the concept of hadronic valence, 
i.e. gg, component of the photon wave function will be introduced \A l25|. 

The main difficulties to be dealt with are: i) how to construct the photon-hadron coupling 
when a gg pair is produced by a photon with g+ > 0; and ii) how to describe the nonvalence 
content relevant for the process under consideration, both in the pion and in the photon 
wave functions. 

The first issue is addressed by using a covariant generalization of the vector meson dom- 
inance approach (see, e.g., {2^) at the level of the photon vertex function (see Ref. |23|]). 

As a matter of fact, it is necessary to construct the Green's function of the interacting gg 
pair in the 1~ channel. For the description of the vector meson vertex functions in the valence 
sector we use the eigenfunctions of the square mass operator proposed in Refs. 

The 

simplified version of the model that we are going to use j|3l| includes confinement through 
a harmonic oscillator potential. The model showed a universal and satisfactory description 
of the experimental values of the masses of both singlet and triplet S-wave mesons and 
the corresponding radial excitations j^^, giving a natural explanation of the almost linear 



3 



relationship between the mass squared of excited states and the radial quantum number 
n js^ Is^. Therefore such a relativistic QCD-inspired model for pseudoscalar and vector 
mesons retains the main feature of the spectra and at the same time allows one to perform 
simple numerical calculations. 

The second issue, i.e. the contribution of the nonvalence (2q2q) components of the pion 
and photon wave functions, is addressed using a model where a quark in the valence com- 
ponent radiates a pair by a contact interaction This interaction is described through a 
pseudoscalar coupling of quark and pion fields, multiplied by a constant. In a recent study 
of meson decay processes within LF dynamics this approximation was shown to give 
a good description of the experimental data. Here, we just follow the above suggestion to 
parameterize the radiative pion emission amplitude from the quark. 

Another important point to be treated carefully is the contribution of the instantaneous 
terms, which is strictly related to the fermionic nature of the constituents. We remind 
the reader th^t the Dirac propagator can be decomposed using the light-front momentum 
components |^], as follows: 



k'^-m^ + te k+{k- - k-^ + ^) 2k+ 

where 7^ = 7° + 7^ and k~^ = (|k_Lp + m'^)/k^. The second term on the right-hand side 
of Eq. Q is an instantaneous term in the light-front time, related to the so-called zero 



modes. As already known (see, e.g., 2l|), the instantaneous contributions play a dominant 
role in the description of the pion electromagnetic form factor in the space-like region, in a 
reference frame where g"^ > 0. Therefore a special care is devoted in the present work to the 
treatment of the instantaneous contributions in the light-cone representation of the fermion 
propagators. In particular the contributions of the zero modes are under control, thanks to 
the momentum behavior of the hadron vertex functions. It should be pointed out that the 
effects of the instantaneous terms is emphasized by the small mass of the pion. 

Our description contains a small set of parameters: the oscillator strength, the constituent 
quark mass, and the width for the vector mesons. We use experimental widths for the 
vector mesons, when available ^|, while for the unknown widths of the radial excitations 
we use a single width as a fitting parameter. The constant involved in the description of 
the nonvalence component can be fixed by the pion charge normalization in the limit of a 
vanishing pion mass. The evaluation of the instantaneous vertex functions involves a further 



parameter (see Sect. X). 

Previously, the elastic time-like form factor was explored in the light-front quantization 
in a boson model of qQ mesons with point-like vertexes jl^ , which does not exploit the rich 
structure of the meson excited states. Here, by studying in a common framework the pion 
space- and time-like form factors, we also access information from the radially excited vector 
meson wave functions. Indeed, in our approach, in the time-like region the virtual photon 
couples directly to the vector meson resonances, which in turn decay in tt+tt". Therefore our 
microscopical model could represent a useful tool to address the investigation of the vector 
meson Green function. 

In the present paper in order to simplify the numerical calculations, we use a massless 
pion, i.e. we evaluate the pion form factor at the chiral limit. In the time-like region the full 
result for the pion form factor is always given by the pair-production process (" Z-diagram" ) 



alone, independe ntly of t 
contribution 



lis approximation. In the space-like region only the "Z-diagram" 
22, 0, 0] survives for a massless pion [2^. The importance of 
the "Z-diagram" contribution to the electromagnetic current for g+ > was also recently 
investigated in the context of the Bethe-Salpeter equation within the light-front quantization 
in Ref. H- 

This work is organized as follows. In Sec. II, we present the general form of the covariant 
electromagnetic form factor of the pion in impulse approximation, which is our starting 
point, and our vector-meson-dominance approach for the dressed photon vertex. In Sec. Ill, 
we first decompose the triangle diagram in on-shell and instantaneous contributions. Then 
we integrate for g+ > on the light-front energy in the momentum loop of the triangle 
diagram, under analytical assumptions for the vertex functions. 

The valence components of the light-front meson and photon wave functions are defined 
in Sec. IV. In Sec. V, we discuss the contribution of the nonvalence component of the photon 
to the time-like current, which appears through the vertex for the radiative emission of pions 
by a virtual quark inside the photon. In this section, we also discuss the contribution of the 
nonvalence component of the pion wave function to the space-like current. In Sec. VI, we 
introduce the pion and vector meson wave functions in the expression of the triangle diagram. 
The pion time-like and space-like form factors written in terms of the valence components 
of the meson wave functions and of the emission/absorption vertexes are derived in Sec. VII 
and Sec. VIII, respectively. In Sec. IX, we briefly revise the light-front model for the pion 
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and vector mesons and conclude the derivation of our model for the pion form factor with a 
discussion on our treatment of the vertex functions for the instantaneous terms. In Sec. X, 
we present the numerical results for the pion form factor in the momentum transfer range 
between -10 {GeV/cY and +10 {GeV/cf'. In Sec. XI, our conclusions are presented. 



II. COVARIANT EM FORM FACTOR OF THE PION 



Our starting point is the covariant expression for the amplitude of the processes tt 7* ^ vr' 

stic case or an antipion in the productior 
(see the triangle diagram of Fig. 1). In 



or 7* — s> tttt', where the meson vr' is a pion in the elastic case or an antipion in the production 
process, evaluated in impulse approximation 'Q 



the time-like region one has (see Fig. 2) 

f = (;r7r|J^(g)|0)= -t2e^N,J-0^A,{k-P^,P,)K{k,P^) x 

Tr[S{k - P^) 7^ S{k - q) r^(fc, q) S{k) 7'] , (3) 

where iVc = 3 is the number of colors; S(p) = is the quark propagator with m 

P — m + le 

the mass of the constituent quark; g'^ is the virtual photon momentum; P^ and P^* are the 



pion momenta. The factor 2 stems from isospin algebra, since 



Tr 



where (1 + r2)/2 is the isospin factor of the current and the other isospin factors in Eq. Q 
pertain to the pions. 

The function AT^{k,Pj^) is the momentum component of the qq vertex function for the 
outgoing pion, which will be taken as a symmetric function of the q, q momenta. In this 
vertex function, Pt^ is the momentum of the outgoing pion and k is the momentum of the 
incoming quark (see Fig. 2). The "bar" notation on the vertex function labels the adjoint 
Bethe-Salpeter amplitude, i.e. the solution of a Bethe-Salpeter equation where the two- 
body irreducible kernel is placed on the right of the amplitude, while for the Bethe-Salpeter 
amplitude it is placed on the left js^, Q • This is a well known property of time orderings 
implied by the Mandelstam formula, for initial and final states (27L l37| . The vertex function 
is defined by the following equation 

f^z _ _|_ i^iz _ jy^z _|_ (27r)^ J 

expz(fc'-2/-A;-a;) (0|T[g(x) g(i/)]|P^) , (5) 



where q{x) is the quark field. 

To obtain the current matrix element for the space-like region, should be replaced by 
— and 7f by n'. Then the pion vertexes A^(A;, P^) and A^{k — P^, P^) in Eq. (jH)) are to be 
changed with A7r(— /CjPtt) and A^^i^k + Pt,,Pt,i), respectively (see Fig. 3). The momentum 
dependence of the vertex functions A^{k — P.„,P^) and A7r(A;, P^) is expected to regularize 
the integrals of Eq. Q. 

The dressed photon- vertex, T'^{k,q), is related to the photon Bethe-Salpeter amplitude, 
which is defined from the three-point function in the standard form: 

1 



k^ — + It fc^ — -|- le 

expi{k' ■ x' — k ■ X + q ■ x")R'^{x, x' , x") . 

The three-point function is given by 



(2vr) 



d X d X d X X 



P^(x,x',x") = (0|T[g(x) q{x")-i^'q{x")q{x')] |0) 



(6) 



(7) 



which is the matrix element between the vacuum states of the time ordered product of the 
four quark fields written above [38]. 

The central assumption of our paper is the microscopical description of the dressed photon 
vertex, r^(A;,g), in the processes where a photon with g+ > decays in a quark- ant iquark 
pair. In these processes we use for the photon vertex, dressed by the interaction between 
the qq pair, the following covariant vector meson dominance approximation (see Fig. 4) 

fvn 



V^{k,q) 
where 



-9^" + 



q^q^ 

n 



Vnv{k,k- q) An{k,q) 



q^-Ml + iMnVniq"] 



(8) 



-9^" + 



qi'q'^'l 



1 



(9) 



,s the vecto. „,eso„ propagate. Q. In » . the deeay eon.ta,. of the „th 

vector meson in a virtual photon, M„ the corresponding mass, A„(A;, q) gives the momentum 
dependence and Vnu{k, k — q) the Dirac structure of the VM vertex function, while r„(g^) 
is the total decay width. If we approximate Eq. (jH)) considering on-shell quantities for the 
VM in the numerator, i.e. if we replace q~ with P~ = (|q_Lp + M^)/g"'', we have 

[^\{Pn)]* fvn 

n,A q'-Ml+tM,,tr,{q'] 



r^(A;,g) = v^5][eA(P„)-V;(A;,A;-P„)] A„(A:,P, 



(10) 
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where the quantity ex{Pn) ■ Vn{k, k — P„) A„(A;, P„) is the VM vertex function and ex{Pn) 
the VM polarization. Note that the total momentum for an on-shell vector meson is = 
{P- = (Iq^p + M^)/q+,Pn± = q±,P+ = while = {q-,q±,q^} and that at the 
production vertex, see Fig. 4, the light-front three-momentum is conserved. 

In Eq. (jHJ the sum runs over all the possible vector mesons. The vector meson decay 
constant, fvn, can be obtained from the definition jl^ 



5^v^/y,„=(O|g(O)7^g(O)|0„,A) 



(11) 



with \(f)n,\) the vector meson state. A detailed expression for fvn is given in Appendix A. 
The total decay width in the denominator of Eqs. and (fTUI) . r„(g^), is vanishing in the 
SL region. In the TL region it is assumed to be equal to 



M2 

n 

g2 



1/2 



(12) 



41 



42|. 



where p(g^) = [g^ — 4m^]^/^/2 

In Ref. [4^ the following expression was used for the Dirac structure, Vn{k, k — Pn), of 
the vector meson vertex : 



vj^ik, k') = r 



ki" + k'l" 



2 Pn-k + mMn 
where k' = k — Pn- 

Let us consider, instead of Eq. (fT^ . a symmetric form for Vn{k,k — Pn) 



(13) 



vj:{k,k-Pn 



k>' + k'>' 



r- 



(14) 



Pn-k- Pn-k' + 2mMn + 2m 

If in Eq. (II 4j) both the CQ's are taken on their mass shell (i.e., k~ = = (|k_Lp + m^)/A;"'") 
and the VM mass, M„, is replaced by the free mass, Mq, of the quarks in a system of total 
momentum g'^. 



Mo = [{kon + (g - k)on) ■ {kon + (g " fc)on)]^^^ 



one obtains 



Vi:{k,k-q) 
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, - (g - kt,^ 

Mn + 2m 



(15) 



(16) 



This form coincides with the on-shell expression given in Ref. |8| for the "^Si vector meson 
vertex, but then T^[k, g) of Eq. ()1U|) is not anymore a four vector. 
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Let us note that to obtain the pion form factor one actually needs only one of the 
components of the current. In the following we will derive the pion form factor from the 
plus component. 

In Ref. 01 to evaluate the pion form factor we considered the plus component of Eq. 
(|lUp. where A„ and the VM polarizations were taken at the vector meson pole, and the 
on-shell expression for V^, as given by Eq. (fT^ . was used, in order to have the structure of 
the VM vertex suggested by the Hamiltonian LF dynamics. 

In Appendix B, starting from Eq. (jH)), we propose a current which satisfies current 
conservation. In the reference frame, where g"*" = M„ > and q_L = 0, the nth term of this 
current has exactly the same plus component as the nth term of the current defined in Eq. 

One might wonder that a bare coupling term should be added to the current defined 
in Eqs. or pO|) . However, as it is shown in Appendix C, in the case of a massless pion a 
bare coupling produces violation of current conservation. Therefore we do not consider this 
term in the present paper. 

III. APPROXIMATING THE TRIANGLE DIAGRAM ON THE LIGHT-FRONT 

Our aim is to retain the essential physics contained in the triangle diagram (Fig. 1) and 
at the same time to construct a bridge toward the Hamiltonian wave function language. 
At the same time we wish to go beyond a simple valence description. To accomplish these 
goals and to eliminate the relative light-front time between the quarks, we perform the k~ 
integration in Eq. Q with some assumptions on the analytical structure of the A and F 
vertexes for the pion and the photon. To be more precise, Eq. (jS)) is evaluated with the 
assumptions that: (i) the momentum components, A{k,P), of the vertex functions, both 
for the pion and the vector mesons, vanish in the complex plane k~ for — * oo; and (ii) 
the contributions of the possible singularities of A(fc, P) can be neglected. Furthermore, the 
Dirac structures of the vector meson vertex function, V^(fc, k — q), are assumed to be regular 
functions of the complex variable k~ . The expressions for VJ^{k,k') given by Eqs. (fTB|) or 
(|Tl|l obviously fulfill the requirement that no pole is present in the A;" complex plane. 

To make clear the discussion of the k^ integration, it is helpful to first separate instanta- 
neous and non-instantaneous contributions, using the decomposition of the Dirac propagator 
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given in Eq. (j2I). Indeed this decomposition is useful to have a better control on possible 
divergences both in k~ and in fc+ integrations. In particular, as already mentioned, it should 
be pointed out the tight relation between the instantaneous terms and the so-called zero 
modes, where = 0. We assume that the behavior of the functions A(k, P) in k^ is able 

Since three propagators are present in Eq. Q, one has a total of eight contributions. The 
contribution with three instantaneous terms vanishes because of the property 7"^ 7+ = 0, 
since the combination 7"*" 7^ 7"*" appears. Also the three contributions with two instantaneous 
terms vanish, as a consequence of our assumptions on A{k,P). Indeed, only a single pole 
from the propagators is present in these contributions. Then, since we assume that the 
functions A{k, P) go to zero for \k~\ -^00 and disregard their singularities, we can perform 
the integration in the k~ complex plane closing the contour in the semiplane where no 
singularity in the propagators is present and we obtain a null result. Moreover, two of these 
contributions with two instantaneous terms are also identically vanishing because of the 
presence of the combination 7"*" 7^ 7"*". 

Therefore we are left with four contributions : three contributions with one instantaneous 
term only and one contribution with no instantaneous term. 

To evaluate the triangle diagram we treat separately the time-like case and the space-like 
case. 



A. Time-like case 



In the time-like case, one has q'^ = P^ + P^, and > . Equation Q written in light- 
front variables becomes (the Jacobian for the transformation to the light-front variables is 
1/2): 

' ~ \2^Yn ' J (k^ - p.-)kHk^ -q-) ^'^^ J 

As-(A; — -Pjr, -P#)A^(A;, P^) /'i7'\ 
(k- - + ^)ik~ -q~-ik- g)-„ + ^)(/c" - P- - (k - P.)-„ + ^^) ^ 

The on-mass- shell values of the minus-components of the momenta in Eq. (fTTj) are given by 
_ ki + m^ _ (k-q)i+m^ (k-P^iW 

%n - ) {'^ Don - ^+ _ ^+ ' ^^)on " ^+ _ p+ ' ^^°) 
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and the operator is defined as follows 



= {f-p^ + m)-f\f-i + m) r^(A:, q) + m) 7^ . (19) 

Let us perform the decomposition of the propagators in instantaneous and in on-shell 
parts (see Eq. 0), as discussed at the beginning of this section. Then Eq. (fTTj) becomes 

3' = JL + Jr + J2+J^ (20) 

where J^^ represents the on-shell contribution and Jtii = 1, 2, 3) represent the contributions 
with one instantaneous term. Then we have 



and 



where 



^ Tr[[(^ - /^^),„ + m] 75[(|^ - ^)„„ + m] r^(A;, g) (|^,„ + m) 7^] 

°" (A:- - A;- + ^)(fc- - - (fc - g)-„ + ^)(fc- - P" - {k - P.)-„ + p^)^ ^ 

^ Tr[7+ 75 [(^ - + m] r/-(fc, q) (|^,„ + m) 7^] 
^ 2(fc--A:- +^)(fc--g--(fc-g)-„ + ^) ^ ^ 

^ Tr[[(^ - /^^)o„. + m] 75 7+ r^(fc, g) (^,„ + m) 7^] 
^ 2{k--k-^ + ^){k--P--{k-P^)-^ + ^) ^ ^ 

^ Trm - p^)on + m] 7^ [{^ - ^),„ + m] T>^{k, q) 7+ 7^] .261 
3 2(fc--g--(fc-g)-„ + ^)(A:--P--(fc-P.)-„ + ^) • ^ ^ 

In Eqs. ()211 the three propagators of the triangle diagram generate three poles: 

, _ _ I.- _ 



k+ - g^ 



7 



If 

Within our assumptions on the vertex functions, A{k,P), if A;+ < there are no poles 
in the lower complex semi-plane of k~ (cf. Eq. ()27|)). Therefore, if the k~ integration is 
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performed by closing the contour of integration in the lower semi-plane, a vanishing result 
is obtained. Furthermore, if > , there are no poles in the upper complex semi-plane 
and a vanishing result is obtained by closing the contour in the upper semi-plane. Then, the 
integrals (j21l \22\ have contributions only for in the range < A;^ < g"*". The integration 
range can be decomposed in two intervals, < A;^ < and < k^ < . In the first 
one, if the k~ integration contour is closed in the lower semi-plane, only the pole /c^-j falls 
within the integration contour, while in the second one, if the integration contour is closed 
in the upper semi-plane, only the pole fc^-, falls within the integration contour. Then in 
the range < A;+ < P+ one has contributions from JTq'^, Ji and J2 , while in the range 
P^ < k^ < one has contributions from J'^^, and . 

Let us introduce the free mass, Mo(fc"*", k_|_; P+, P_l), of a qq system of mass M, total 
momentum P, kinematical momenta (A;"'", k^) and {P^ — k^, P± — k^) for the pair of quarks: 



Mo2(fc+,k^;P+,Pj 



k]_ + (P - k)5_ + 



(28) 



X 1 — X 

where x = k^/P^, with < a; < 1. Using this definition of the free mass, the following 
equations hold: 

1 P+ 



(P- _ (p _ k)-^ 



k- 



(M2 + P^ 



P^ 



(M2-M2(A;+,k^;P+,P^)) 



(29) 



1 



P 



[p/- _ (p/ _ (fc _ p))-^ _ (A; - P)-J + p. 

p/+ 



(P-k)2^+m2 (P'-(k-P))2^+m2 



(30) 



[M2 - M2((A; - P)+, (k - P)^; P'+, P'^)] 
where x' = {k+ - P+)/P'+. 

Then, performing the k~ integration and using Eqs. from Eq. we obtain 



e m" 



dk~^dh 



{2nr Jo {k+-P+)k+{q+-k 

The quantities If and I2 in Eq. (jSH) are defined as follows 



{0(p+ - k+) If + e(fc+ - p+) (31) 



Ii — A-^[k, P„)Ajf[k P^, P^^ 



.^on,(l) + ^1,(1) +^2,(1) 



(32) 



-'2 



A7r(A;, Pn-) A^(A; Pj^^P^^ 



-{k-q)o„ 



^on,{2) -^1,(2) -^3,1 



(2). 



(33) 
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where 



rptJ. 

on,(l) 



Tr[[(^ - p^)on + m] 75[(^ - ^)„„, + m] r'^(l) (Iton + m) 7^ 



[M2(fc+,k^;g+,qx)-g2 



mt 



on,{2) 



+ m] r'^(2) (|£o„ + m) 7^] 



[M2(fc+, k^; g+, q^) - - te] ml - Mi{{k+ - P+), (k - P^)^; P^+, P^j 



rpfl 



Tr[7+ 75[(|^ - + m] r^(^) (^,^ + m) 7' 
2 [M2(A;+,k^;g+,q^)-g2_,e] 



1,2) 



rpH 
^2,(1) 



P^ 



^3,(2) 



PJ 



Tr[[(|^ - /^.),„ + m] 75 7+ r/-(l) + m) 7^] 
2 [M2(A;+,k^;P+,P^j 



+ m] 1/^(2) 7+ 7^] 



Mo2((^+-^+),(k-P.)±;P^+,P^±) -m2 



with 



r^«=r^(fc+,k^,fc" = A;^),g) (2 = 1,2) 



(34) 
.(35) 

(36) 
(37) 
(38) 



(39) 



The first term of Eq. (jHT|) . with fc"*" — P+ < 0, and the second term, with fc"*" — P+ > 0, 
are represented in Fig. 2 by the diagrams (a), with the arrow of k~^ — from the left to 
the right, and (b), with the arrow of k'^ — P^ from the right to the left, respectively. In 
the first term only the vertex function A^(fc, P^) has the momentum fraction k'^ /P^ in the 
valence- sector range [0, 1] and in the second term only the vertex function Ajf{k — P^^, Pjt) 
has the momentum fraction {k^ — P^)/P^ in the valence-sector Tange [0, 1]. 



B. Space-like case 

In the space-like case, one has P^, = P^ + q'^. Then the expression for the triangle 
diagram can be obtained from Eq. Q replacing — P^ with P^, tx with tt' and the pion 
vertexes A^(fc, P^) and A^(A; — P^^, P^) with !s.T,{—k, P^) and k.„i{k + P^, P-k')-, respectively : 



m 



-i2e—N, 



d^k 



Tr[Sik + P^h'S{k-q) V'ik^q) S{k)-i'']K'{k + P^, P^,)K{-k, P^) 



/2 V (27r)4 

^ dk~dk'^dk 



e m 

TT-rrrA^, 



(27r)4/r-^ J {k+ + P^)k+ik+-q+) ^ J 

A,,(fc + P^,P^OA.(-A;,P,) 



{k- - k-, + ^){k- -q--ik- g)„-, + pp^)(A;- + P- - (A: + P,)" + 



(40) 
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where 



O'^ = {^ + p^ + m)^\^-i + m) T^{k, q) + m) 7^ (41) 

As in the time-hke case, let us decompose the propagators in on-shell and instantaneous 
parts. Then Eq. (pn|) becomes 

f = J'L + J"i + J"2 + J", (42) 

where J'^ represents the on-shell contribution and J'^{i = 1,2,3) represent the contribu- 
tions with one instantaneous term. Then we have 



and 



^ ' (27r)4 fl'^^ J [k+ + P+)k+{k+-q+) ^Vlft + /^.,/^.'JiWl ft, ^^.J . I44j 

where 

^ rr[[(^ + + m] 75[(^ - ^)„„ + m] r^(A:, g) + m) 7^] 

°" (A:- - A;- + ^)(A:- - g- - (A: - g)„- + ^){k- + P- - (A; + P.)-„ + ^ 

^ ^ Tr [7+ 75 [(1^ - ^),„ + m] r^(fc, q) {^^ + m) 7^] 
' ' 2(A:--fc- +^)(fc--g--(fc-g)- +^) ^ ^ 

^ Tr + p^),^ + m] 7^ 7+ r^(A;, g) + m) 7^] 
2 2(A;--fc-„ + ^)(A;-+P--(A: + P0-„ + ^) ^ ^ 



^„ ^ Tr m + A)on + m] 7^ [(1^ - i)^ + m] r/-(fc, q) 7+ 7^] . . 

3 2{k--q--{k-q)-^ + ^){k-+P--{k + P^)-^ + ^) ■ ^ > 

In Eq. (jlOj) the quark propagators generate three poles : 

ft(l) — fton ^+ 5 

K2) = 1' + {k-q)on - 



^("4) = + + ^^)on - , (49) 
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where 



{k + P.) 



(k 



(50) 



k+ + P+ 

Let us assume > 0. Therefore, if fc^ > and within the hypotheses stated at the 
beginning of this Section, there are no poles in the upper complex semi-plane and a vanishing 
result is obtained by closing the k~ integration contour in this semi-plane. Furthermore, if 
k^ < —P^i there are no poles in the lower complex semi-plane of k~ . Therefore, one obtains 
a vanishing result by closing the contour of integration in the lower semi-plane. Then, the 
integral has contributions only for A;+ in the range —P^ < k^ < . The integration range 
can be decomposed in two intervals, —P^ < k^ < and < < In the first one, 
if the integration contour is closed in the lower semi-plane, only the pole fc^-j falls within 
the integration contour, while in the second one, if the integration contour is closed in 
the upper semi-plane, only the pole k'^2) f^^s within the integration contour. Then in the 
range —P^ < k'^ < one has contributions from and J'"^, while in the range 

< A;+ < g+ one has contributions from JT'^^, J'"^ and JT"'; 
As a consequence, can be decomposed as follows 



(51) 



where j^^^^ has the integration on k^ constrained by —P^ ^ k^ < 0, while j^^^^f^ has 
the integration on k~^ in the interval < /c^ < g^. As illustrated below, the valence 
component of the pion contributes to j^^^^ only, while 

is the contribution of the pair 
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22 



23 



production mechanism from an incoming virtual photon with g+ > ^ , 3, 
35l |. Performing the k~ integration, the two contributions to are given by the following 
expressions 



J 



m 



(2^)^ 



dk~^d]<. I 



P+ {k+ + P+) k+ {q+ - k+) 



^on,(4) + ^2,{i) + -'3,(4) 



X 



-^■K/{k + PjT) PlT/)-^1T{~k, Pjt) 



k — ^P-K +(^+-P7r)o7i 



(52) 



e m^,^ /■'?+ 



dk~^d]i 



(27r)3 /2 Jo {k+ + P+) k+ (g+ - k+) 



^on,{2) + -'1,(2) + 



(2) 



A7r/(fc + Pyr, P7r/)A7r( k,P.^^ 



k-=q-+{k-q)on 



(53) 
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where 



rpl^l 

on, (4) 



[P- - + + P.^ - (A: + + (A; - g). 



(54) 



^/^ _ Tr [[(1^ + + m] 75[(^ - + m] r^(2) + m) 7^] 



■ on,{2) 



[q- + {k- g)-„ - A;-„ + ze] P^, - (fc + P^)-„ + (A; - g)- 



(55) 



_ Tr [7+ 75[(|^ - 4),^ + m] r/-(2) (^,„ + m) 7^] 
''^'^ 2 [q- + {k- g)-„ - fc- + Z6] 



(56) 



rri/fl 

^2,(4) 



Tr [[(^ + f^),„ + m] 7' 7+ r^(4) (|^on + m) 7 
2 [P- - (A; + P.)-„ + /.-J 



(57) 



^3,(4) 



[[(1^ + A)on + m] - i)on + m] r^(4) 7+ 7^1 



P^, - {k + P^)-^ + {k-q)„^ 



(58) 



with 



^ Tr [[(1^ + /^.).„ + m] 7'[(|^ - ^),„ + m] r^(2) 7+ 7^ 

Ptt/ — (/C + P-K)on + {k — q)or 



r^(4) = r'^(A;+,k^,A;- = fc(-),g) 



(59) 



(60) 



The contributions and j^-^-'^)'^ are represented by diagrams (a) and (b) of Fig. 3, 
respectively. 

1. Valence region contribution 



Let us change integration variables in Eq. ()52|) for the valence contribution, defining 
k'~^ = k~^ + P^ and = k_L + Pttx, with (A;'^, k'j_) the light-front momentum of a quark in 
the valence range. Then j(^)'^ acquaints the following more familiar expression 

^2 rP+ dk'+dk!, 



3 



e m 



(27r)3 /2 " 7o (A;'+ - P+)k'+{P^, - /e'^ 



rjiffl I T""'/^ I n~'^ {-^ 

^on,(4) + ^2,(4) + ^3,(4) 



X 



A,r/(A; , Pn-/)A^(P^ A; , Ptt 



k'- — k'~ ' 



(61) 
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where 



k'+ 



(62) 



and we have defined k + = k' . The quantities T'J^^^-^, Tg^'^^-j, ^3^"^^-^ can now be 
expressed as follows : 

_ r P+ P+ [(^on + m) - + m] r^(4) [{f - p.U + m] 7^] 



^2,(4) 



Tr [(^on + m) 75 7+ r^(4) [(^ - yP.)on + m] 7^] 

2 K-M2(A;'+,kl;P+,P.^)] 



(64) 



rp/fl 

^3,(4) 



-PJ 



Tr [(^on + m) 75[(r - /^^,)o„ + m] T>^iA) 7+ 7^] 



(65) 



In Eq. ()6ip both the vertex functions have the quark momentum fractions k'~^/P^, and 
(Ptt — k')^/P^ in the valence sector [0,1]. Note that the on-shell momenta in Eq. ()63p 
allow one to retrieve the relativistic spin coupling factors with the spin 1/2 Melosh rotations 
automatically included ^, l^l • 



2. Pair-production contribution 



By making use of Eq. the quantities T'J^ f^^) ; '^1%) ' '^3%) pair-production 

contribution (Eq. (jSSl)) become 

Tr [[{^ + p^)on + m]75[(|£ - ^)„, + m] r^(2) {^on + m) 7^] 



on,{2) ^TTf 



q- - qo + le ml - M^{k'+, k^; P+, P^^^^ 



(66) 



Tr [7+ 75[(^ - ^),„ + m] r^(2) (^,„ + m) 7^] 

^1,(2)- 



(67) 



T^'M _ p+ 
^3,(2) — ^TTl 



Tr [[{f + /',)„„ + m] 75[(|^ - i)on + m] r'^(2) 7+ 7^] 



m 



TT/ M^i^k'^ , k^ ; P^, , Ptt/x ) 



(68) 



where go = ^on + ~ ^)on fc'"*" = A;+ + P+, k'_i_ = k_L + P,r±- In Eq- dSSj) only the vertex 
function A,r/ has the quark momentum fraction (P,r + k)^ /Pi in the range [0, 1]. 
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IV. VALENCE COMPONENT OF THE LIGHT-FRONT WAVE FUNCTION 



A. Meson wave functions 



1. Pion 

To fully interpret the terms that appear in Eq. ()3ip and in Eqs. ()52|) . ()53p . we have 
to discuss valence and nonvalence components of the light-front wave function. Let us first 
)egin with the valence component, which can be derived from the Bethe-Salpeter amplitude 
4q . In the present model the pion Bethe-Salpeter amplitude is given by 

m jt + m , ^-p. + m 

where the pion vertex is j- 'y^ATj{k, P^). 

The valence component of the light-front wave function can be obtained from the Bethe- 
Salpeter amplitude (j^^ in the valence sector, < < P+, disregarding the instantaneous 
terms in Eq. (jU^ . multiplying "^t^ by the factor [k^ {k^ — P^)]/{2m) and integrating over 
k- : 

Two poles, k{i) and A;(3), appear in Eq. (jTU)), respectively in the lower and in the upper k~ 
semiplanes. We perform the k~ integration in the lower complex semi-plane disregarding the 
contributions that arise from the singularities of the vertex A7r(fc, Ptt) (cf. the assumptions 
(i) and (ii) at the beginning of Sect. III). Then the pion wave function becomes 

Uk\ k.; Pi, P,J ,n) |(^ - p,U + -1 (71) 

If the k~ integration is done in the upper semi-plane within the same assumptions, one has: 
0^(fc+,k^;P+,P^^) = 

In principle, the elimination of the relative light-front time between the quark and the 
antiquark in the pion Bethe-Salpeter amplitude by the k~ integration in Eq. (|7n|) . should 
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give a unique answer, which defines the valence component of the wave function in the range 
< < P^, with both the quarks on their mass shelL Therefore in order to require con- 
sistency within our model, we will assume [A^(A;, P7r)][fe-=fc-j and [A^(/c, Pn-)][^,_^p-_(.p^_^-|-^-| 
to be equal in that kinematical range (note that M^{k+, k±; P+, P^^) is equal to M^{P+ - 
k^, P^_|_ — k_|_; P+, P,r±))- This assumption produces a momentum component of the valence 
light-front wave function symmetrical for the exchange of the quark momenta, since the 
vertex function A^(A;,P7r) is assumed to be symmetrical. 

Within a Bethe-Salpeter approach, the function (p^^ fulfills a two-body Schroedinger- 

like equation, with the proper Melosh structure represented by the matrix {^on + 

n 

m) 7^ [{jk — pT^)on + '^l [8]. Therefore, when the plus component of the quark momentum 



is in the interval < A;"'' < P^, will be identified in our approach with t 
wave function, with momentum component ?/'7r(fc"'", k_L; P^, P^rx) (see Ref. |8 

0^(A;+, k^; P+, P^^) = (|^,„ + m) 7' U - A)on + m] ^^{k+, k^; P+, P, 



le HLFD pion 
a) . (73) 



2. Vector meson 



In analogy with the pion case, one can define the light-front VM wave function, which 
describes the valence component of the meson state \nX). Indeed, starting from the Bethe- 
Salpeter amphtude for a vector meson 



'^nxik,Pn 



^ + m 



^2 _ ^2 _j_ 



exiPn)-Vnik,k-P^) An{k,Pn 



(74) 



{k — PnY — rn? + le 

the valence component of the light-front wave function can be defined from '^n\{ki Pn) 
integrating over k~ , disregarding the instantaneous terms and multiplying by the factor 
[k~^ {k~^ — P^)]/{2m), as we already did for the pion. Furthermore, in this case one has to 
take on their mass shell both the quark momenta in the Dirac structure of the VM vertex 
function, Vn{k,k — Pn) : 



0nA(A:+,^±;P+,P„j 



-tk+{k^ 



p' 



27r k'^ — im? + le 



An{k,Pn) ex{Pn)-[Vn{k,k-Pr,)i 



+ m 



{k Pn 



(75) 



where \Vn{k,k — Pn)\on is defined by Eq. (fTBj) in order to retrieve the ^Si vector meson 
vertex of Ref. ^. Assuming that An{k,Pn) does not diverge in the complex plane k~ for 
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\k I ^ oo, and neglecting the contributions of its singularities in the k integration, the 
valence VM wave function is 



0„A(A:+,kx;P+,P„. 

'e^{R,)-[Vr,ik,k-P^)] 



In analogy with the pion case, we assume [A„(A;, P„)]j^_^^-^] = [A„(A;, P„)][^_^p-_(^p^_^^-^] in 
the valence sector, < A;^ < P^. 

As for the pion, the function 0„a, with the plus component of the quark momentum in 
the interval < A;'*' < P^, will be identified with the HLFD vector meson wave function, 
with momentum component ilJn{k~^ ,\<.±; Pn ,Pn±) , j^h 

</.„A(fc+,kx;P+,P„a) = 

= i^on + m) [ex{Pn)-[Vnik,k-P^)]on\ ^ - A)on + m] t/-, ( A;+ , k^ ; P+, P,^) .(77) 

In conclusion, Eqs. (ffT| I7!?|) and f7^ I77jl establish a link between the momentum part of 
the meson HLFD wave functions and the momentum part of the meson vertex functions. 

The valence component of the VM wave function are normalized to the probability of the 
valence component of the meson state \nX) (see Appendix D). This probability is estimated 
in a schematic model in Appendix E. 

The corresponding normalization for the pion wave function is included in an overall 
normalization constant for the pion form factor. 

B. Photon wave function 

One can define as well the valence component of the hadronic contribution to the photon 
wave function, starting from the Bethe-Salpeter amplitude of the photon, which can be 
written as: 



^ pi^r.(M)(^^^,^t^ . (78) 



where T^^{k^q) is the photon vertex amplitude (see Eq. (jH))). 

In analogy with Eq. (f73j) . the valence component of the virtual photon light-front wave 
function can be obtained from the Bethe-Salpeter amplitude (fTSj) in the valence sector. 
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< A;"*" < g"*", separating out the instantaneous terms of Eq. (f75jl. integrating in k~ and 
multiplying by the factor [k^ (A;+ — P^)]/{2m). 

Then, using our explicit expression for r^{k,q) given by Eq. (jSI), the light-front wave 
function of the photon can be defined by 

0^;(fc+,k^;g^g^q^) = 

= -^kHk^ - q^) f ^ J- \1 [rnA:,g)],„ 'f^ , (79) 

J 2tt k^ — + te [k — qY — + le 

where the label "on" in [r^(fc,g)]Q„ means that, as in the VM case, the Dirac structures of 
the photon vertex amplitude, r^(fc,g), have to be taken with both the quark momenta on 
their mass shell. Therefore, in analogy with Eq. (f7S|) . a possible choice for the quantity 
Vnik^k — q) of Eq. (jHl) is given by the quantity Vnik^ k — q) as defined by Eq. (fT^ . 

L \ on 

Then, performing the k~ integration with the assumptions given at the beginning of 
Section 3, in the range < /c"*' < Eq. (fTUj) becomes 



0^;(fc+,kx;g2,g+,qx) = (|^o„ + m) 7/'j;:(A;+,kx;g2,g+,qx) [(|^-^)o„ + m] , (80) 

where the function ilj^{k~^, k±] g^, g"*", q^), which includes the Dirac structures of the photon 
vertex, is defined by 

As in the previous meson case, to have consistency in our virtual photon wave function 
model, we will not distinguish between [A„(A;, g)]j^_^^-^] and [An{k,q)]^i^^^^-_(^^_f^^-^-^ in the 
valence sector, < A;"'' < g^. Therefore we obtain for il)^{k^ , k_L; g^, g+, q^) the same result 
when the k~ integration is performed both in the lower or in the upper k~ complex semiplane. 

The valence wave function, 0f^(fc"^, k_L; g^, g"*", q^), and the function ?/'j^(A;+, k_L; g^, g"*", q^) 
depend on the value of g^ carried by the virtual photon. Note that in the time-like case a 
singularity appears in the photon valence wave function (see Eq. ()81|) ). 

If, as in Eq. fllO|) . the photon vertex \r^{k^q)\on is taken with on-shell quantities for the 
vector mesons in the numerator, i.e. if we take 

[r^(fc,g)]„„ = v^E^A(Pn) • P„)l„„ [K{KPn)\,-=,,. r ^5^1^"!". ..-( 82) 



n,\ 



q^-Ml + iMnVn{q'] 
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and identify Eqs. (fTHjl and (fTTjl . then the function ipf^^ik'^ , k_L; q^, g"*", q^) can be written as 
follows: 



X 



n,A 




g2-M2 + ^M„f„(g2) 



(83) 



in terms of the momentum part of the HLFD vector meson wave functions, 
V'n(A;+,k^;P+,P„^). 

V. CONTRIBUTION OF NONVALENCE COMPONENTS TO THE CURRENT 
A. Time-like case: the photon nonvalence component 

The process of pion-antipion production is shown in Fig. 2, where the dashed lines (both 
in (a) and in (b)) represent two different light-front times. At the first time (the one on 
the right) the hadronic valence component of the virtual photon is represented, while at the 
second one the 2g2g photon nonvalence component is depicted (see also Fig. 5). The two 
parts of Fig. 2, i.e. (a) and (b), differ by the emission vertex of an antipion or of a pion (see 
also Fig. 5 (b)), respectively. The corresponding quark amplitudes for the radiation of an 
antipion or a pion are given in Eq. by the antipion vertex K^{k — P^; P^), evaluated at 
k~ = /c~„ for (/c— P^)+ < 0, and by the pion vertex A^(A;; P^), evaluated at k~ = q~ + {k—q)~j^ 
for fc+ > P+, respectively. 

Once the interaction that couples the valence to the 2g2g component is known (see Fig. 
2), the amplitude for the photon decay in a ttW pair can be constructed. To this end, let us 
introduce a kernel operator /C which realizes this coupling. Then, we can write the following 
equation to relate the valence component of the pion wave function, ilj.^{k'^ ,'k'j_; ,P.^±), 
to the vertex function Ajrlk; P^^) at k~ = q~ + {k — g)~„, which is the amplitude for the pion 
emission (see Fig. 2 (b) and Fig. 6): 




(7')^"(7')/.'a' JC-J {k'+,k'^ ; fc+,k^ ; g-,g+,qx) 



(84) 
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For simplicity, the example of a 7^ structure was used in Eq. ()84p . just to be consistent with 
our assumption of a pseudoscalar pion model. 

One can write an analogous expression for the emission of yf : 

ITl 

J TV 

1 ^ fP^ dk'+dk' 



{l'f%l')p'c.' /C^;^' {k'^, kl ; fc+ - P+, - P.^ ; g', g+, q^) . (85) 

In our model calculation, both pion emission vertexes will be substituted by a constant, 
following Ref. H- 



B. Space-like case: the pion nonvalence component 

A part the direct photon coupling to the quark line, in the space-like region the nonvalence 
component of the final pion wave function appears for g+ > in both the contributions of 
the current obtained after the k~ integration, and given by Eqs. (j52p and (j53|) (see Fig. 
7). On one hand the valence component of the final pion is coupled to the nonvalence 
2g2g component (see Fig. 7 (b)), through an interaction kernel 7i, which contributes to the 
quark-photon absorption vertex of Eq. ()52|) . given by r'^(4) = T^{k' — P^^, q) with k'~ = k'^. 
On the other hand the vertex A.„{—k; P^) in Eq. (fS^. evaluated at k^ = + {k ~ q)^^ for 
—k~^ < 0, describes the quark-pion absorption through another interaction kernel JC' and 
generates the nonvalence 2g2g component of the final pion (see Fig. 7 (c)). We identify the 
kernel JC' with the kernel /C, already used in the previous subsection A for the description 
of the pion emission (Fig. 6). 

Equation ()52j) gives a contribution to the SL form factor where the initial and the final 
pion valence components appear (diagram (a) of Fig. 3). The plus component of the quark- 
photon absorption vertex, given by r'^{k' — P^^^q) with /c'" = k'^^ which enters in Eq. (|52j). 
is represented by an empty circle in diagram (a) of Fig. 3 and is approximated by the sum 
of i) the bare photon vertex multiplied by a renormalization constant, a, (diagram (a) of 
Fig. 7) and ii) the contribution due to the 2g2g component of the final pion wave function, 
which is represented by diagram (b) of Fig. 7. 
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Therefore, we can make the following identification: 



^,7o k"+{q+-k"+) 



a'l3' 

As already discussed in Sec. II, we do not consider the bare term photon vertex in the 
present paper, since it violates current conservation for a massless pion (see Appendix C). 
Therefore, disregarding the bare photon vertex in the right-hand side of Eq. ()8fjj) . we can 
formally write: 

[r+(A:'-P.,g)](,_,-^)^H^+. (87) 

One could try to interpret Eq. (jSTjl in terms of constituent quark form factors. However, we 
have to point out that the absorption vertex of Eq. (jHHjl does not depend only on g^, as one 
could naively think, but it depends on the virtuality of the quark, and therefore depends on 
the hadron where this process occurs. 

Let us note that, within our assumption of a vanishing pion mass, the contribution of 
Eq. is also vanishing (see Sect. VIII) and therefore there is no contribution from 

[r+(A:'-P.,g)](,_,-,). 

Equation represents the pair-production term (Z-diagram) and is depicted in Fig. 7 
(c). The quark-pion absorption vertex, given by A.„{—k; P^) evaluated at k~ = OT + ik — q)~^, 
which appears in Eq. (jSHjl can be written as 

il')^aiiy'lC"Jf,, k^; kl ; P^, P+, P.,^) ^k'^, k^; P+, P^^) . (88) 

For our purpose this quark-pion absorption vertex will be taken consta,nt, as we do in the 
TL case for the quark-pion emission vertex, as was proposed in Ref. [l8| . 

VI. TRIANGLE DIAGRAM AND PION LF WAVE FUNCTION 
A. Time-like case 

Let us insert into Eq. ()31|) the photon vertex of Eq. pO|) . Furthermore, whenever the full 
expression for the light-front pion wave function 0^(A;+, k_|_; P+, Pjr±), given by Eq. (fTHl. 
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appears in Eq. ()3ip and the momentum fraction is in the valence-sector range [0, 1], let us 
replace it with the expression of Eq. (fTSj) . i.e. let us write the pion vertex in terms of the 
momentum component of the HLFD pion wave function. This means that we introduce in 
Eq. (jSH) the wave functions ^^{k+, k^; P+, P^^) and V^f ((^^ - ^/), (k - P^)±; P^, P#±) 
when these wave functions have the correct support. Then the triangle diagram can be 
expressed as follows: 



e m ,^ ri 



dk~^d]<.\ 



E 



V2 [e^(P„)]7yn 



{0(P+ - A;+) + Q{k+ - P+) . 



The quantities /i,n,A and /2,n,A in Eq. (j89|) are defined as follows 



h,n,X = [k^{k- P^,P^) kn{k,Pn\.^^- 

[ [ 

Ug2-M2(A;+,k^;g+,q^)+ze] 



X 



(89) 



TonXl,n,\) + ri,(l,n,A) + r2,(l,n,A) \ (90) 



/2,n,A = [K{k,P^) An{k,Pn) 



k =q +(k-q)orL 



X 



[g2-M2(A;+,k^;g+,q^)+ze] 



Ton,{2,n,\) + 7'i^(2,„,A) + ^3,(2,n,A) \ (91) 



where 



Ton,{l,n,A)= C(^+,k^;^.+ ,P^±) X 



(92) 



Tr 



[(^ - A)on + 75[(^ - ^)„„ + m] eA(P„) ■ Vn{K k-Pr 



T, 



on,{2,n,\) 



V^;((fc+-P+),(k-P.)^;P+,P^^) X 



(93) 



Tr 



[(^ - A)on + m] 7^[(^ - ^)o„ + m] eA(P„) ■ V;(A;, k - P„) 



fc =g +{k-q)„. 



_ {^on + m) 7^ 



2"l,(l,n,A) 



1 m 

2 U 



Tr 



7^ 7^[(l^-^)o„ + m] exiPn)-Vn{k,k- Pn) (^on + m) 7 



(94) 
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1 in 

Tl,(2,n,X) ^ ~ 2 Y ^^^^^ ~ -^'^'-^*)]fc-=g-+{fc-g)-„ 



Tr 



7+ 7'[(^-^)on + m] ex{Pn)-Vn{k,k-Pn 



k-=q- + (k-q)o„ 



(95) 



^2,(l,n 



(l,n,A) 



Tr 



[(1^ - A)on + m] 7^ 7+ eA(P„)-K(A;,A;-P, 



(^on + m) 7^ 



(96) 



1 



■3,(2,n,A) 



C((^+-P+),(k-P.)x;P^,P. 



X 



(97) 



Tr 



[(1^ - P.)on + m] 7^[(|^ - i)on + m] k(Pn) " %{k, k - P„) 



7 7^ 



k =q +{k~q)on 

Let us notice that the momentum component of the LF pion wave function does not appear 
in the instantaneous terms Ti (i „ a) and Ti (2,„,a)5 because in these terms the propagator 
[{f, — pT,)on + m]/ {k~ — — {k — Pt,)^^ + 1+^%+ ) is replaced by 7"^/2. Indeed the amphtude 

1 m 



2 U 



7 



(98) 



does not obey the same two-body Schroedinger-hke equation as the hght-front pion wave 
function (j)T,{k'^, k±; P+, Pttx) does. 

As aheady noted at the end of Sect. Ill A, the first and the second term of Eq. ()89j] are 
represented in Fig. 2 by the diagrams (a) and (b), respectively. Note that, due to the 9 
functions, the final pion or antipion wave functions enter into the first or the second term of 
Eq. (jHUj) . respectively. In Eq. (jH^ the pion vertexes [Ajf{k — Pt,] Pw)], evaluated at k~ = k~^, 
and [A7r(fc; Ptt)], evaluated at k~ = q~ + [k — q)~n, have the momentum fraction outside the 
valence-sector range [0, 1] and can be related to the quark amplitudes for radiative antipion 
or pion emission, respectively (see Figs. 2 (a), 2 (b)). The presence of these vertexes gives 
rise to the contribution of the nonvalence component of the virtual-photon wave function, 
relevant for the process under consideration. In the space-like region the analogous processes 
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can be interpreted for g"*" > as the contribution of the nonvalence component of the pion 
wave function in the final state These points have aheady been illustrated in Section 
V. 

If we choose to take the Dirac structures in the photon vertex with both quarks on their 
mass shell, i.e. T^{k^q) = [r^{k,q)]^^ (see Eq. (j82j) ). then whenever the full expression 
for the light-front vector meson wave function 4>n\{k^ , k±; , Pn±), given by Eq. (fTUj). 
appears in Eq. (jH!?j) . we can take advantage of our identification of Eqs. (fTBj) and (f77j) 
to express the quantities /i,n,A and l2,n,x through the momentum component of the HLFD 
VM wave function. However, in the instantaneous contributions to /i,„,a and l2,n,x which 
are proportional to the quantities T2 (i „ a) and T3 (2^„^a) "we do not express the VM vertex 
functions [A„(A;, P„)]^_^^- and [A„(/c, P„)]^_^^_^j^^_^^- through the momentum component 
of the HLFD VM wave function, because the full expression for this function given by Eq. 
(I77|) does not appear in these instantaneous terms. 

Then we obtain 

f Mk^, k^; ^n+, Pn±) [Ml - M^{k+, k^; P+, P„^)] 



(l,n,A) 



+ 



\ [q^-M^{k+,kr,q^,q±)+^e] 

[An{k,Pn) ],_,-^ T2,(i,„,A)} (99) 



'2,ra,A 



A^(/c, Pn-^ 



k-=q- + {k-g)o„ 



X 



kr, P+, Pn±) [Ml - Ml{k+, kx; P+, P„^)] 
[g2-M2(A;+,k^;g+,q^)+ze] 



,n,A) 



+ 



[An{k,Pn)] 



k =q +(k-q) 



o„ ^3,(2,n,A)} 



(100) 



The quantities To„^(i,„,A), Ti^^i^^^X), ^2,(i,n,A) and the quantities To„^(2,n,A), ^i,(2,n,A), ^3,(2,n,A) 
in Eq. and in Eq. have the same expressions as in Eqs. (jH^ . (jMj) . (PHjl 

and in Eqs. (PHj). ((HHI), and (P7|). respectively, with eA(Pri) ■ t4(^, ^ — -Pn) ^_ ^_ and 
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exiPn) ■ Vnik, k - P„) _ both replaced by ex{Pn) ■ V^ik, k - P„) (see Eq. 

-i K —q -ryK—q)on -Ion 

(ITH|) for the definition of \VJk, A; - P„) 1 ). 

L J on 

Note that the region of integration over k'^ in Eq. ()89|1 (a consequence of the non vanish- 
ing integration in the k~ complex plane) agrees with the support of the wave functions ipT^ 
and if)^ of Eqs. ()92|) and (jUHj) . respectively. Furthermore, in agreement with the above as- 
sumptions, the vertex associated with the virtual photon and consequently the wave function 
ipn{k^, k±; P^, P„_l) in Eqs. and (jl(JU|) have the intrinsic fraction of the plus-component 
of the quark momentum, k'^ /q^ = k^/P^, in the interval [0, 1]. 

To be able to evaluate the TL pion form factor we have still to assign a value, in the instan- 
taneous terms, to the VM vertex functions [A„(A;, P„) and [An{k, Pn)]i^-^^-_^_^i^_^y , 



as well as to the pion vertex functions A^(A;; P^ 



and [A^f(A;- P^,P^)]; 



=<? +(k-q)on' 



B. Space-like case 

As in the time-like case, let us replace in Eqs. (jUT|) and ()53p the pion vertex function 
with its expression in terms of the momentum component of the HLFD pion wave function, 
whenever the full expression for the LP pion wave function, (Eq. (jZ^) appears, taking 
advantage of our identification of Eqs. ()71|1 and (ffH|) . 



1. Valence region contribution 

Substituting in Eq. (jfilj) the pion initial and final wave functions and noting that for the 
final pion the "bar" vertex gives the complex conjugate wave function, while in the initial 
state the vertex gives the initial pion wave function, in the valence region one gets (see Fig. 
3): 



(27r) 



dk'+dk\ 



{k'+-P+)k'+{P+-k'^ 



+ Tj(4) z^^(P+-A;'+,P^^-kl;P+,P^^) A^,(A;',P^,) 



+ 
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+ T 



3,(4) 



(101) 



where 



r 



on, (4) 



Tr 



ilt'on + m) ^'li^' - p^,)on + m] r^(4) W - r.)on + m] 7' 



(102) 



T: 



2,(4) 



Tr 

2 



^'on + m) 7' 7+ r^(4) [(r - r.)on + m] 7^ 



(103) 



^ 3,(4) 



- Tr 
2 



+ m) 75[(^' - r.^,^ + m] r'^(4) 7+ 7^ 



(104) 



In the instantaneous terms proportional to and to T3 (-4-) we do not express the pion 



3,(4) 



vertex functions 



,_ and [A^(P7r — fc', ^77)]^/-=^'- in terms of the momentum 



component of the HLFD pion wave function, because of the presence of 7+ instead of [(^' — 
Ar/)on + "i] or [(1^' - f^)on + , respectively. 

In Eq. (imri the photon vertex, r''(4) = V^{k' - Pn,q), evaluated at k'- = k'^, is the 
amplitude for the photon absorption by a quark. As discussed in Sect. V B, the photon 
absorption operator can be decomposed in a bare vertex, i.e. 7^^ (Fig. 7 (a)), plus other 
terms. From the expansion in the light-front Fock-space, the next term relevant for the 
process we are analyzing is due to a contribution of the nonvalence 2g2g component of the 
final pion wave function, see diagram (b) of Fig. 7. This contribution can be thought of as 
an expectation value of an operator between the valence component of the wave functions 
for the initial and final pions. The operator can be constructed by applying to the virtual 
photon wave-function the kernel, 7i, which produces the nonvalence pion component from 
the valence one (see Eq. 



2. Pair-production contribution 



Also the pair-production contribution to the current, can be rewritten in terms of the 
momentum component of the HLFD pion wave function (see Eqs. ()71|1 . (ffH|l ) when the 
light-front pion wave function appears. Then Eq. becomes (see Fig. 3 (b)): 

■HI), ^ _eN^^ K c/fc+rfk^ [A.(-A:;P.)],_^^_^(,,^)--^ ^ [e^(Pn)]*/yn 

^ (27r)3/Jo (fc+ + P+) A;+ (g+ - A;+) 

( n+ 

[A„(A;, Pn)\k-=q- + {k-q)- 



q 



[q^-Miik+,kr,q^,q±)+^e] 



^on,(2,n) + ^l',(2,n) + ^3,(2,n) f (105) 
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where 



Tr 



U + A)on + "^]7^ [(1^ - i)on + m] eA(P„) " Vn{k, k - P„) 



-{k-q)on 



(106) 



{2,n) 



1 m 

2 A 



k =q +{k-q)on 



(107) 



Tr 



k =q +(k-q)on 



T' 



C,((fc+ + ^.^),(k + P.)±;Pj;,P./±) X 



(108) 



Tr 



U + P.)on + m] 7^ 



+ m] 



ex{Pn)-Vn{k,k-Pn] 



7^ 7^ 



fc =q +{k-q)on 

As for the time-hke case, we have used the expression of Eq. |T0|) for the photon vertex 
with the virtual photon going into a qq pair. The "bar" vertex A^r/ imphes that the final 
pion wave function in the above expressions has to be complex conjugated. If we take the 
Dirac structures in the photon vertex with both the quarks on their mass shell, as in the 
time-like case (see Eq. (jH21)), then using Eqs. (fTBj) and (f77j) we can express through 
the momentum component of the HLFD vector meson wave functions, when the LP VM 
wave function is present, i.e. in the terms given by Eqs. and (jlOTj) : 



k =q +{k-q)o„ 



E 



V2 [6^(P„)]* /. 



Vn 



{2nfhJo {k+ + P^) k+ {q+ - k+) [g2-M2 + zM„f„(g2 

>„(A;+, k^; P+, P„x) [Ml - M^ik+, k^; P+, P„^)] 



[g2-M2(fc+,k^;g+,q^) + 



^on,{2,n) + ^l,(2,ra 



(2,n) 



+ 



[An{k, Pn)]fc-=q-+(fc_g)-^ ^3,(2,n)} 



with 



(109) 



on,{2,n) 



rA{k+ + P+),{k+P^)r,P:-nPn,±) X 



(110) 



Tr 



[(^ + A)on + m]7^ [(1^ - 4)on + m] ex{Pn) ■ Vn{k, k - Pn) ilton + m) 7 
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rpl 



1 m 



2 h 



k =q +(k-q)o„ 



X 



Tr 



7+ 7' [(1^ - i)on + m] exiPn) ■ Vnik, k - Pn) 



+ m)7^ (111) 



TUn) = \ rUik^ + P;^), (k + P.)x; P+, P.,±) X (112) 



-3,(2,n) ~ 2 

Tr 



[(1^ + A)on + m] 7^ [(1^ - 4) on + m] ex{Pn) ■ Vn{k, k - Pn) 



1^ 7' 

on 



The vertex A^(— A;; P^) evaluated at /c"" = g~ + (A; — g)~„ represents the pion absorption 
amphtude by an on-shell quark. The presence of this process can be also interpreted as a 
2q2q component in the final pion wave function (see Figs. 3 (b) and 7 (c)), as illustrated in 
Sect. V. 

VII. TIME-LIKE EM FORM FACTOR OF THE PION 

We have pointed out in the Introduction that, for a unified description of TL and SL 
form factors, it is necessary to choose a reference frame where the plus component of the 
momentum transfer, q"*", is different from zero (otherwise, = q~^q^—(l\ cannot be positive). 



Therefore, as in Ref. l26||, in order to calculate the pion form factor we adopt a reference 
frame where q_L = and g"*" > 0. 

The decay of a time-like virtual photon is written in terms of the time-like form factor of 
the pion as follows 

r = (7r7r|g(0)7^g(0)|0) = e {P^^ - P^) F^{q') (113) 

where = P^ + P^ is the four momentum of the virtual photon. In Fig. 2, the diagram- 
matic analysis of the virtual-photon decay in a tttt pair is shown. 

The virtual-photon decay amplitude can be obtained from Eq. pi3|l by evaluating the 
plus-component of the matrix element, j^. To be able to evaluate the matrix element 
from Eq. (|H^. we substitute in Eq. (jH^ constant values for the vertexes, and V^, 
namely for pion or antipion radiation by a quark, Eqs. (j8l|l and (j85|l . respectively (see also 
Eqs. (jniH IHT jl ). Then, it remains to specify the values of the instantaneous vertex functions 
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inEq. jni, [A^{k-P^,P^)l^ 



in Eq. JHSI), and [A„(fc,P„)] in 



k-=k-„ ' ' " ' ■" =1 +(k~l)on 

Eqs. (jnnmnni)- TWs win be thoroughly discussed in Sect. IX. 

By using Eqs. and (jHn|l one can obtain the pion form factor FT^{q^) from the 

plus-component of the current : 

fvn 



(114) 



[g2-M2 + 2M„f„(g2) 

where 5'y„(i?^), for > 0, is the form factor for the VM decay in a pair of pions, as expected 
from the VMD approximation. The characteristic feature of our approach is that we aim at 
a microscopic description of both /y„ and gvniQ"^)- 

Let us now evaluate J2x of Eq. (jHI?|) . The momentum of the vector meson is = {P~ = 
(|q±P + M!^)/q~^, q_L, g+} and the momentum of the virtual photon is = {q~ , q±, q~^} (as 
already noted, see Fig. 4, at the production vertex the LP three-momentum is conserved). 
In a frame where = 0, one has q~ = q'^/q^ for the photon, while for the vector meson 
P~ = M'^/q'^. Using Cartesian components for the four vectors, i.e. = [a°,a], in this 
frame the three polarization four-vectors are given by 



e^^ = [0, 1, 0, 0], = [0, 0, 1, 0], = [P„./M„, 0, 0, ^1 + r,] (115) 

where r] = P^^/M^. . Let us recall that, in the frame we are adopting, Pnz = {q'^ — P~)/2 = 
(g"*"^ — M^)/2g"'". Therefore, in the reference frame defined by q_L = and > the 
polarization four- vector does not have a defined sign for the zero component. The plus- 
component of is given by 

' _2 



Mr. 



+ 



\ 



1 + 



P2 

nz 



-Mi 



+ 



1 + 



g^ 



M2 



2g+M„ 



(116) 
0) 



2g+M„ 

The plus- component of the other polarization four- vectors are vanishing (i.e., = 
and therefore we have Ea 4 (^n) ^\{.Pn) ■ % = [et{Pn)T (^z{Pn) ■ %■ 

Each term of in Eq. ()114|) is invariant under LP boosts, that are kinematical, and 
therefore to simplify the calculations it can be evaluated in the rest frame of the correspond- 
ing resonance. In the rest frame of the nth— resonance one has g"*" = M„ and g~ = g^/M„ 
for the photon, while P^ = P~ = Mn for the vector meson. This means that we choose a 
different frame for each resonance, but all these frames are related by kinematical LP boosts 
along the z axis to each other, and to the Breit frame where g"*" = — g" = v^— g^, adopted in 



previous analyses of the SL region (we have always qj_ = 0) |2lL |26| 
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Then, in the evaluation of the sum in Eq. ()114p one has = 1 and 

Y.[tUPn)\ tx{Pn)-Vn = -Vnz 
A 

for the contribution of any resonance! In conclusion we have 



(117) 



1 



Nr 



9+ V2 dk+d\^A 



{e(P+ - A;+) + e(A:+ - P+) h^J^ , (118) 



10 {k+ - P+)k+{q+ -k\ 
where the quantities Ji,„ and l2,n can be obtained from Eqs. (|^. ()1UU|) replacing 



m//^ [A^(A;-P^, P^f) ]^_^^-^ and m//^ [A^(A;,P^) 
tively (see Eqs. and 



k~=q +{k-q)o 



with and P^, respec- 



• k^; P+, 0) [Ml - M^jk^, kr, P+, 0] 

[g2-M2(A;+,k^;g+,0)+ze] 



%n,(l,n) + '^,(l,n) 



+ 



[A„(/c,P„)]fc_^fc-^ T2,(i,„)} 



(119) 



■2,n 



1 [q^-MUk+,kr,q+,0)+te] 

[^n{k,Pn)]k-=q-+(k-g)-„ 7^,(2,™)} 



Z5n,(2,n) + ^,(2,ra) 



+ 



(120) 



with 



on,(l,n) 



- V^:(A;+,k^;P+,P.^) X (121) 
Tr [[(1^ - p^)on + m] -f^[{^ - i)on + m] \Vn,{k, k-P^)] {^^n + m) 7' 



on,{2,n) 



C((fc+-P.+),(k-P.)x;P/,P.±) 



X 



(122) 



Tr 



[(1^ - Ar)on + 7^[(|^ - ^)on + Vnz{k,k - Pn) {^n + m) 'J 



Tr 



7"" 7'[(^-^)on + m] Vn,{k,k-Pn) (^0^ + ^)7^ (123) 
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1 in 



X 



Tr 



r l%f-i)on + m] Vnz{k,k-Pn) (^on + m) 7' 



(124) 



2,(l,n) 



X 



Tr 



[(1^ - A)on + m] 7^ 7+ Vnz{k,k-Pn) {hn + m) 



(125) 



p.) , ; P 
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7r±J 



Tr 



+ m 



Vnz{k, k Pn) 



7^ 7' 



(126) 



The ^5*1 vector meson vertex [t4(fc, k — Pn)]on given by Eq. (fT^ . as it was used in previous 
calculations 8], is completely determined by the kinematical momenta of the individual 
quark and antiquark. In the ^5*1 vector meson intrinsic frame, where q'^ = Mn and = 0, 
one has 

1,3 _ ( p3 _ t.3\ \ 
'^on l-* n Jon \ 



Vnz{k, k Pn) 



f - 



Mon + 2m 



(127) 



where M^^ = (jk^p + m^)/x{l - x) {x = k+/q^ = k+/Mn). 

The kinematics for the final two-pion state in the particular frame where the photon has 
momentum = 0, = M„ and = q^ /Mn can be derived from the energy-momentum 
conservation, which yields 



and thus 



+ P: 



P+ + p^ 



P + ml 



(128) 



P + ml 



1 IP, 



P + ml 



P^ 



(129) 



Pjt g+ (1 - x^) 

where = P/^ /q'^ and a;^ = P^ jq^ = 1 — a;,^- Eqs. (I128|) and (jl29p put in evidence the 
relation between the kinematical variables of the virtual photon and the ones of both pion 
and antipion. In the time-like region the value of does not fully determine the values for 
the four-momenta of the pion and the antipion in the final state of the tttt pair. In order to 
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reduce the freedom we make the purely longitudinal choice, i.e. P^f_L = — PttX = 0. Then, 
from Eqs. ()128|) and p29|) . one obtains 

Let us note that the minimum allowed value for is 4m^. At this threshold value one has 
= 1/2. Therefore, since 

- P+ = g+ - 2P+ = g+(l - 2x^) = M„(l - 2xJ , (131) 

one cannot evaluate Eq. ()118j) exactly at threshold, unless an exact cancellation occurs 
between vanishing numerator and denominator of Eq. ()118j) . For finite values of m^r, the 
values = 1 oi = Q are possible only for an infinite value of the momentum transfer 
and imply an infinite value of P^^ or P^^i respectively. 

In the limit of a vanishing pion mass {tjIt^ = 0), Eq. p30|) gives x^^ = 1 or 0, which 
implies that one of the terms of Eq. pi8j) vanishes due to the function. To simplify 
our calculations, in the following we make the approximation m^r = and adopt the choice 
Xt, = 0, which implies = 0, = 07= l^n, = = Mn, and P^ = 0. Then only 
the second term of Eq. pi8|) . containing the quantity l2,n, gives a contribution to the TL 
pion form factor. 

Furthermore, for = one has '7^n,{2,n) = . Indeed, in this limit — pn)on = ^on 
and 

{fon + m) 7^ {fon + m) = {fon + m) (- fon + m) 7^ = (- fonhn + m^) 7^ = (132) 

Therefore only the instantaneous contributions 7j,(2,n) and ^,(2,n) survive in the limit of a 
vanishing pion mass and can be written as follows : 



1 in 

Tr [7+ lilt - i)on + m] Iv^k, k - P„) 



X 



+ m) 



(133) 



7/>;(A;+,k^;M„,0j 



X 



Tr 



[-hn + m] [{f: - (^)on + m] Vnz{k, k - P, 



(134) 
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To evaluate the time-like pion form factor we have still to specify the values 
of the instantaneous vertex functions [Aff(fc — Ptt, -P7r)]fe-=g-+(fc_g)- in Eq- P33|) and 
[A„(A;, P„)]^_^^__i_|.^_^^- in Eq. (|12(jp that, as already explained, cannot be directly related 
to and ipn- 



VIII. SPACE-LIKE EM FORM FACTOR OF THE PION 

The space-like form factor of the pion can be obtained from the plus component of the 
proper current matrix element 

f = {7r\mi'q{0)H) = e {P^ + P^,) F^{q') (135) 

where = P^, — P^ . 

In our reference frame, where q± = and > 0, the minus- component of the four- 
momentum transfer is given by q~ = q'^/q^, which is negative in the space-like region. Let 
us note that 

(136) 

Hence, the constraint < is obviously fulfilled for any value of Pn±, since |P7r/±| = \Ptt±\ 
and P+ = g+ + P+ > P+. From Eq. (IT!?H|l one has 

' " ^ P.^ + P.^) " (1 + X.) ' ^ ^ 

where = P^/q^. Therefore, once a value for |P7r±| is chosen, P^ and Pj^ are fixed. For 
a purely longitudinal motion of the pions, i.e. Pn± = Pit/± = 0, it is easy to obtain from 
Eq. (Unil) that 

In the limit ofrrin = 0, the longitudinal momenta of the pions according to Eq. p38j) are 

P+ = and P+ = g+ (139) 

for any value of the momentum transfer. 

In a frame where q~^ ^ 0, the electromagnetic current in the space-like region, Eq. 
(|3T|). receives contributions from the valence component of the wave function, j*-^-*"*' given 
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by Eq. (jlUip . as well as from the nonvalence components, j'^^^^^ of Eq. (jlU9|) . i.e. from the 
Z-diagram contribution (see Fig. 7). 

The contribution of the pion valence wave function to the current can be calculated 
from Eq. p01|) introducing the plus component of the operator T^{k — Ptt^q) ior k~ = 
as given in Eq. (jHSI) and discussed in Sect. V B, once the values of the pion vertex 



functions 



and [A7r(P7r — A;', P7r)]fc'-=fc'- the instantaneous terms have 



been specified. 

In the limit of zero pion mass, according to Eq. (jl(Jl|) the valence contribution to the 
space-like pion form factor vanishes, since = 0. Then, only the Z-diagram contribution 
survives in this limit, as in the time-like region. 

The contribution of the Z-diagram to the elastic pion form factor can be obtained from 
Eq. ()1H5|1 by substituting in Eq. the pion absorption vertex of Eq. ()88|1 . The result 

can be written as follows: 

fvn 



(140) 



q' - Ml 

Since f^iq^) is invariant under kinematical LF boosts, we choose to evaluate the contribu- 
tion of each vector meson, fn{q^)-, in the same reference frame that we used in the time-like 
region, i.e., we adopt the rest frame for each resonance (g"*" = M„, q_L = ; = g+ = M„, 
P- = Ml/q+ = M„). 

Then for a finite value of the pion mass we have: 

et dk+ 



fn\q') 



with 



T' 



7r3 P+ + P+ Jo k+ {q+ - k+) (P+ + k+) 



' V;„(A:+, k^; P+, 0^) [Ml - Mjjk^ , k^; P+, 0^)] 
[g2-M2(A;+,k^;g+,0^)+^e] 

[^"(^'^")]fc-=g-+{fc-g)o„ ^3,(2,n)} 



CX(A:+ + P+),(k + P.K;Pj;,P.,±) X 



T' -I- T' 

■'■on,{1,n) ^ l,{2,n) 



+ 



;i4i) 
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Tr 



[{f + A)on + m] 7^ [(^ - i)on + m] Vnz{k, k - P„) {^on + m) 7 



1 m 



{2,n) 



2 U 



A^/(A; + Pn-, Pj^ 



=q +(fc-c?)o„ 



X 



Tr 



7+7M(^-^)on + m] VUk,k-Pn) 



+ m) 7^ 



(143) 
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T' 

3,(2,n) 



V^:,((^+ + ^.+),(k + P.K;P+,P 



X 



Tr 



[(1^ + A)on + m] 7^ [(1^ - + m] V;^(fc, A; - P, 



7 7 



(144) 



The Dirac structure, Vnz(k, k — Pn) , for the ^Si meson state is given by Eq. (11271) . As 

L J on 

already noted, in our reference frame one has = 1. Equations p4H|l and ()144j) represent 
the instantaneous contributions. Analogously to the time-like case, for a vanishing pion 

mass one has T^^f^^^n) = (^^^ -^l- 

Let us now evaluate f^^iq"^) at — * 0^ for a finite value of the mass of the pion. To 
begin witlLj_we consider: i) a constant value for V^^, ii) a simple form for the LF pion wave 
function [201 



V^.,[(A;+ + P+),(k^ + P.^);P+,P. 



m 



PJ 



and iii) in the instantaneous term p4H|l ) take A^,(A; + P^,P^ 



7r_L 5 TT/ 5 7r/_L 



:(145) 



proportional 



fc =9 +(fc-(j)o.. 

to ipl,[{k^ + P^), (k_L + Pn±)', P^i, Ptt/x] (see the next Section). For a finite value of the 
mass of the pion, let us note that in the limit — 0~ from Eq. p38j) one obtains oo 
and P^, = {Mn + P^) oo. Then, since the squared free mass for the final pion is 



Mol,P+ + P+),(kj 



) . p+ p 



PJ 



+ 



|k_Lp + 



,(146) 



^ P+ + k+ p+-p+-k+^ 

it becomes large for P+ ^ oo, i.e. M^, ~ P+. Then ^:,[(A;+ + P+), (k^ + P^^); P+, P^,^] 
becomes a constant for P^ — oo. Furthermore, for P/ — oo the traces in Eqs. 
are proportional to ~ P^. Therefore, collecting together the factors P^ in Eq. 



one 



concludes that for a finite value of the pion mass lim 



/i^(g2) ~ lim 



1/P+ = 0. 



The same result, lim ^2^0- fi (q ) = 0, should also hold for pion wave functions which are 
eigenf unctions of a Hamiltonian |31|. 

On the contrary, in the limit of = 0, the longitudinal momenta of the pions are 
P/ = and Pj; = Mn, respectively (see Eq. fll38p ). Then, according to Eq. p01|) . the 
valence contribution to the space-like pion form factor vanishes, while the Z-diagram yields 
a nonzero contribution. 

A comment is appropriate here. In the work of Ref. j2l|, where 7^ 0, it was found 
that the wave function contribution to the space-like pion form factor strongly decreases 
in the frame = \/ —q^ as the momentum transfer — increases. As a consequence, 
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the Z-diagram contribution, which is zero at g2 = 0, becomes the dominant one at high 
momentum transfer. As the pion mass is artificially decreased in that model, we find that 
the momentum at which the Z-diagram starts to dominate the form factor tends toward 
zero, in agreement with the previous discussion. 

Since in this paper we work at the chiral limit of a vanishing pion mass, in our reference 
frame, the full space-like pion form factor is given by F^\q^) alone. It has to be noted that, 
as occurs in the time-like region and for the same reasons, for = only the instantaneous 
terms ^'(2,n) ^^"^ '^3,(2,n) -^is- ^I143p and p44|) ) give contribution to the pion form factor. 
These terms can be written in the following form : 



T' 



1 m 

2 U 



k =q +{k~q)on 



X 



Tr 



7^ [(^ - i)on + m] Vnzik, k - Pn) 



+ m) 



(147) 



T' 

3,(2,n) 



1 



Tr 



[-fton + m] [{It - ^)on + m] Vnz{k, k - P, 



:i48) 



IX. A LIGHT-FRONT MODEL 

To evaluate the pion form factor we need : 

i) a model for the HLFD pion and vector meson wave functions which appear in Eqs. 
(fTTH|l and (HUD; 

ii) a value for the probability, Pqq,n, of the VM valence component (see Appendices D 
and E); 

iii) an approximation for the pion vertex functions which represent the pion emission or 
absorption by a quark; 

iv) to assign a value to the pion and VM vertex functions with an instantaneous quark 
leg. 
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The vector-meson resonances are described by an effective light-front model inspired by 
QCD jsil, that can be also applied to the pion. The squared mass-operator for the S'-mesons 
contains a harmonic oscillator interaction featuring the confinement and a Dirac delta- 
function that acts in the ^5*0 channel (with a renormalized strength). The wave functions 
for the ^5*1 channel are solutions of the following eigenvalue problem 



4(|/€p + m^) + -^tuV^ + a 



34 J *r(r)=M^vl/^«(r) , (149) 

where = Mq/4 — is the square of the intrinsic quark three- momentum, = 
n uj + Mp and the eigenfunctions \l/^*^(r) are the three-dimensional harmonic oscillator 
wave functions for zero angular momentum. The HLFD wave functions, without the Melosh 
rotations, are given in the Fourier space by 

Mk^, k^, P^, Pn±) = ""(l/^n • (150) 

The factor comes from the different normalizations used for '?/'„(A;"'", k_L, P^, P„_l) and 
\l/^'^(|/tp). Indeed the function \E'^'^(|/tp) is normalized through the equation 

^^''{\K\')\'d^K = l , (151) 



while the function ipn{k~^ ,'k±, ,Pn±) is normalized through the evaluation of the charge 
form factor of a vector meson at = 0, i.e. by using the so-called charge normalization (Ap- 



pendix D), more appropriate in a relativistic context j22|. In the actual calculation, we have 
to consider that, after properly integrating the valence component, its probability should be 
recovered. This amounts to construct a schematic model for the probability, Pn,qq, for each 
excited state ( see Appendix E), and subsequently to renormalize ?/'„(fc"'", k_L, P+, P„_l) in 
Eq. (fTHn|l as follows 



Mk^,^±,Pn,Pn±) = JPn,q,Prl^n\>^\') (152) 



In the model of Ref. j31| the complete form of the pion wave function is an eigenstate of 
the mass operator of Eq. ()149|) plus a Dirac-delta interaction (in the configuration space), 
which is necessary for producing a pion with a small mass (i.e. a collapsing qq pair in the ^5*0 
channel). The pion wave function is found from the pole of the resolvent, explicitly written 
in Ref. j^^]- The result is the following: 

ml 
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where \I/^'^(0) is the S-wave harmonic oscillator eigenfunction in coordinate space at the 
origin. 

In this model, the pion wave function approaches the asymptotic limit, Eq. (jl45|) . imposed 
by the presence of the Dirac delta-function in the interaction. 



The relativistic constituent quark model of Ref. [31[ achieves a satisfactory description 
of the experimental masses for both singlet and triplet S'-wave mesons, with a natural 
explanation of the "lachello-Anisovitch law" js^ Is^, namely the almost linear relation 
between the square mass of the excited states and the radial quantum number n. Since the 
model does not include the mixing between isoscalar and isovector mesons, in this paper we 
include only the contributions of the isovector p-like vector mesons. 

As already discussed in Sec. VI, we approximate the pion vertex functions which represent 
the antipion and pion emission by a quark, as well as the quark-pion absorption vertex by 
means of a constant 

V^ = V^= — X^- and V^ = — X^ (154) 

in agreement with the constant form proposed in Ref. and successfully tested in the 
study of the pseudo-scalar meson decays. The actual value of the constant A,r is fixed by 
the pion charge normalization. 

As anticipated in Sec. VII and Sec. VIII, to simplify our calculations we are going 
to use = 0. Within this assumption, for the time-like form factor only the instan- 
taneous contributions 7^,(2,n) and ^,(2,n) survive, while for the space-like form factor only 
the instantaneous terms '7^'(2n) ^^"^ '^'(2n) remain. Then to fully evaluate the pion form 
factor in the time-like and in the space-like region we have still to assign a value to the 
pion and VM vertex functions with an instantaneous quark leg, i.e. to the vertex functions 
[K^{k-P^,P^\. and X,(A; + P,,P,,)1 in Eqs. (USSl) and (HHj), 

respectively, and to the vertex function [A„(/c, P„)]j^._^g__|_(.^_^-j-^ of Eqs. ()134p and ()148p . 
The instantaneous contributions to the time-like pion form factor corresponding to the ver- 
tex functions [A^(A; - P^, P^)]^_^g__^(^_g)-^ and [A„(/i;, P„)]^_^g__^(^_g)-^ are represented by 
diagrams (a) and (b) of Fig. 8, respectively. 

Let us note that the presence of the factors {k'^ ^Pn) and in the denominators of the 
two instantaneous terms produces an enhancement of the contributions around the values 
{k~^ ± P^) = and /c+ = in the k~^ integration. Within our assumption of a vanishing 
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pion mass, this means that, for both the instantaneous terms, there is an enhancement of 
the contribution at the end point = 0, which corresponds to an infinite value of the 
z component of the intrinsic quark three-momentum, = Mq{2x — l)/2 {k^ = ~oo for 
X = 0, since Mq —>■ oo). Therefore the high momentum part of the meson vertex functions, 
i.e. the short-range part in coordinate space, is very relevant. Then in the vertex functions 
with an instantaneous quark leg, A^*^*^-), we assume that the very short-range part of the 
one-gluon-exchange interaction, which includes spin-spin terms 47], is the dominant one. 
In symbolic notation we have (see Fig. 8) : 

where /C*^* is the Bethe-Salpeter kernel for the instantaneous vertex function A'^*, Gq the 
propagator of two free quarks and A-^"'' the full vertex function. 

The kernel /C*** is assumed to be dominated by the short-range part of the interaction. 
Actually we drastically simplify Eq. (|155|) as follows : 

This amounts to naively assume that A-^"" is an eigenstate of /C*''* Gq. Furthermore, we 
assume that A-^"" is still related to the LF meson wave function as illustrated in Sec. IV, 
i.e. Af^') = [M^(„) - M2]/P4). 

The constant c is thought to roughly describe the effects of the short-range interaction. 



In particular, i 
results of Ref. 



' we take grossly into account only the spin-spin interaction term, then the 
23[ are recovered i) by choosing c = —3/4 for the pion vertex function (Fig. 



8 (a)) and c = 1/4 for the VM vertex function (Fig. 8 (b)) and ii) by using the probabilities 



1 



Pqq-n = — ^ /^^ 3 fo^ the VM valeuce components with 6 uj2 /2 = 1 (see Appendix E). With 
this choice for the constants c's, the relative weight of the VM instantaneous terms with 
respect to the pion instantaneous terms is equal to —1/3. At variance with Ref. in the 
present paper we use this relative weight, wvm = cym/ctt, as a free parameter. 

In conclusion, we replace the momentum component of the pion vertex function in Eq. 
as follows 
Tn 

^^(A;+ - P+, - P,^; P+, P^^) [ml - M^{k+ - P+, - P^^; P+, P^x)](157) 



P 
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and in Eq. ()147p as follows 



m 



k-=q +(k-q)on 



-^rUk^ + P;^, + P^±; P;^n P^/±) [ml - Ml{k+ + P+, + P^x; P+, P.,x)](158) 



The momentum component of the VM vertex function in Eqs. ()1H4|1 . ()148|l is approximated 

by 



CVM 



V>„(A;+, k^; P+, P„^) [Ml - M^{k+, k^; P+, P„x)^159) 



[A„(/c, P„)]^_^g__^(^_^)-^ ^ 

As explained in the previous sections, in the limit of a vanishing pion mass both in the 
time-like and in the space-like case one has P^ = and Pjj = P/ = M„. Then, the 
quantities gyn{.<f) of Eq. ()118|) and f^i^f) of Eq- f)141|) acquire the same functional form, 
despite the sign of g^, and reduce to the same function ini'f') '■ 

dk+ 



^-^^ ) 16vr3 Ml Jo {k+Y (M„ - k+) 



dk I 



Ti,„(A;+,kx) + r3,„(A;+,kj 



X 



<,(fc^,kx;^n,Ox) M2-M2(fc+,kx;M„,0x) ^„(fc+, k^, M„, 



where 7i^„ and Ts „ are given by 



^ ^ _ K-Mo^(fc+,k^;M„,0^)] 
[g2-M2(fc+,kx;M„,0x) + ie] 
Tr 7+ [(1^ - + m] K^l^;, A; - Pn) {l^on + 



(160) 



- 4 



K-M2(A;+,k^;M„,0x)] 



[g2_M2(A;+,k^;M„,0x)+ze] 

k^{k — q)on,z -\- {k — (l)on ' kon + (^^ ~ ^n)kon,z ~ ^1 

- m (2A;+ - Mn){kon - {q - k)on)z Hs{M^) 



(161) 



= WM Tr [[-|^on + m] [{f - (^)on + m] \ynz{k, k - Pn) \ 7 ^ 

= WvM 4 [-A;+(A; - g)^^,^ + {k - q)on ■ kon + {k'^ - Mn)kon,z - m? 

+ mMr, [kon -{q- k)on], Hs{Mo)] . (162) 
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In the last steps in Eqs. fll61|) and p62|) the traces have been exphcitly evaluated and the 
function Hs{Mq) is given by: 

Hs(Mo) = • (163) 

^ ^ Mo + 2m ^ ^ 

Actually the value of c^r together with the value of A^r is fixed by the charge normalization 
and we have to assign a value only to the relative weight wvm- 

In Eq. fll60|) there is no divergence from the poles at the end points = and — k^ = 
0, because of the Gaussian decrease of the VM wave functions at these end points, which 
correspond to infinite values of the z component of the intrinsic quark three-momentum, 
Kz = Mq{2x — l)/2 = — oo or Kz = +oo for a; = or a; = 1, respectively). 

Finally, both in the time-like and in the space-like regions, the pion electromagnetic form 
factor can be written as 



We stress that the pion form factor is continuous at = in the limit m^r — >■ and that only 
the instantaneous terms contribute in this limit. We would like to remind the reader that 
the vector meson wave functions are normalized to the probability of the valence component, 
which can be roughly estimated in a simple model, as shown in Appendix E. The decreasing 
probability of the valence component for the excited vector meson states is essential to make 
convergent the sum over the resonances. 



X. RESULTS 

The pion electromagnetic form factor is calculated through Eqs. ()l(i4|l and pfiOj) . where 
the pion and vector meson wave functions are eigenstates of the square mass operator defined 
in Eq. ()149p ( shown for the vector channel only). 

In our calculation we have a small set of parameters: i) the constituent quark mass, ii) 
the oscillator strength u, iii) the widths for the vector mesons, r„, and iv) the relative weight 
Wvm of the two instantaneous contributions. 

The up-down quark mass is fixed at 0.265 GeV j^^l and the oscillator strength is fixed 
at = 1.556 GeV^ of Ref. |3l|. 

For the first four vector mesons, the masses and widths, presented in Table I, are used. 
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The non-trivial dependence of ^niff) in our microscopical model allows a small shift 
of the VM masses with respect to the values obtained in the analyses of the experimental 
data by using Breit-Wigner functions with constant values for ^n{(f)- 

For the radial excitations with M„ > 2.150 GeV , the mass values corresponding to the 
model of Ref. ^ are used. For the unknown widths we use a single width as a fitting 
parameter. We choose the value r„ = 0.15 GeV , which presents the best agreement with 
the compilation of the experimental data of Ref. 5^. We consider 20 resonances in our 
calculation to obtain stability of the results up to = 10 {GeV/c)'^. 

The probabilities Pqq,n of the valence component of the VM states are fixed according to 
the schematic model of Appendix E (see Eq. ()E16|) and Table II). 

As we discussed in the previous Sections, it is also necessary to know the amplitude for 
the virtual process where a constituent quark radiates or absorbs a pion. This unknown 
function was first investigated in a phenomenological study of decay processes within LF 
dy„a„.i. Q. it was app.oxi„,ated by a constant, obtaining a satisfactory desc.pt.on of the 
experimental data. We followed the approximation proposed in [1^ in the calculation of 
the decay amplitude ^n(Q'^) of Eq. (jlfiOj) . The value of the constant A^r, together with the 
constant 0,^ (see the previous Section), is fixed by the charge normalization. 

The values of the coupling constants, fvn, are calculated using Eq. ()A4|) of Appendix 
A from the model VM wave functions. The corresponding partial decay width, Fg+e-, for 
these mesons are calculated from our values of fvn using Eq. ()A5|1 |lfil | and are reported 
in Table I. The partial decay widths for the vector mesons are in good agreement with the 
data, when available Q]. 

We perform two sets of calculations, to test the effect of the pion wave function model. 
In one set we use the asymptotic form of the pion valence wave function, Eq. ()145|) . and in 
another one we choose the eigenstate of the square mass operator of the model of Ref. j^|, 
given by the pion wave function of Eq. ()153|) . 

The results for the form factor are shown in Figs. 9, 10, and 11. In Fig. 9 the results cor- 
responding to the weight wvm = —0.7 are shown, while in Fig. 11 the results corresponding 
to Wvm = —0.7 and wvm = -1-5 are compared in a linear scale around the p meson peak. 
.n Fig. 9 we also report the results calculated with the masses and the widths used in Ref. 
23[ and reported in Table III. For this calculations the oscillator strength uj = 1.39 GeV"^, 
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the probabilities Pqq-n = —7^ — f and = —3/4, cvm = 1/4 have been used. 

Let us note that our results are the same within a few percent, if in Dirac structure of 
the nth VM vertex (Eq. the free mass is replaced with M„. 

In Fig. 9, we show our results in a wide region of square momentum transfers, from -10 
up to 10 (GeV/c)\ comparing them with the data collected by Baldini et al. 5^ and with 
the data of Ref. 53[. A general qualitative agreement with the data is seen in this wide 
range of momentum transfers, independently of the detailed form of the pion wave function. 
It has to be stressed that the heights of the TL bumps directly depend on the calculated 
values of fvn and ^niQ^)- 

The results obtained with the asymptotic pion wave function and the full model present 
some difference only above 3 {GeV/c)"^. 

The pion form factor is particularly very well described in the space-like region, both 
using the weight wvm = —0.7 or the weight wvm = —1-5, as can be clearly seen in Fig. 
10, where the ratio of the SL form factor to the monopole factor M(g^) = 1/(1 — g^/M^) 
is shown. The excellent agreement with the experimental form factor at low momentum 
transfers is expected, since we have built-in the generalized p-meson dominance. 

The time-like region between and 3 (GeV/c)'^, where p(770), p(1450) and p(1700) 
appear, is shown in Fig. 11 in a linear scale. The p-meson peak is placed at the right 
position using a bare mass of 770 MeV. From this figure it is clear that the parameter 
Wvm is able to control the region of the p(770) peak, while in other regions its effect is less 
relevant. For wvm = —1.5, the p(770) peak is very well described, except for the region 
around 2 {GeV/c)'^, where our results underestimate the experimental data. This dip is 
due to a destructive interference between the contributions of p(770), p(1450), and p(1700), 
and could be potentially sensitive for a detailed test of the model presently adopted for the 
meson wave functions and other approximations introduced. 

It is clear that the introduction of cu-like and 0-like mesons could improve the description 
of the data in the TL region. However, a consistent dynamical description of the mixing of 
isospin states is far beyond the present work, and we leave it for future developments of the 
model. 

Finally, we have also calculate the adimensional quantity, Ft^^q"^) g^/.77^, up to = 
— 1000 (GeV/c)'^, observing a smooth decreases from a value of 0.691 for = —100 {GeV/cY 
to a value of 0.677 for = -1000 {GeV/cy. 
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XI. SUMMARY AND CONCLUSIONS 



In this work, we are able to give a unified description of tlie pion electromagnetic form 
factor in the space- and time-like regions, thanks to the choice of a reference frame where 
g+ > 0. 

The main steps are shortly summarized. Within the covariant approach proposed by 
Mandelstam 2^, the matrix elements of the electromagnetic current are evaluated between 
pion states, in impulse approximation, but with all the vertexes of the triangle diagram 
properly dressed. Exploiting a suitable decomposition of the fermionic propagators, one 
singles out on-shell and instantaneous contributions. The integration over the light-front 
energy, k~ , in the momentum loop of the triangle diagram is performed disregarding the 
effect of possible singularities of the vertex functions and taking care of only the singularities 
in the propagators. 

For the photon vertex function, in the processes where a gg-pair in the odd-parity spin-1 
channel is produced, we use a generalization of the Vector Meson Dominance approach, built 
up from the VM Bethe-Salpeter amplitude (phenomenologically determined) and the VM 
propagator, enlightening the relation between the hadronic part of the photon valence wave 
function and the pion electromagnetic form factor. 

The obtained expression for the electromagnetic current matrix elements are carefully 
discussed and the different contributions are interpreted in terms of valence and nonvalence 
components of the pion and photon wave functions. 

In the valence components of the pion and VM amplitudes, the momentum part is de- 
scribed through the corresponding HLFD wave functions, evaluated in a QCD-inspired model 
which shows a satisfactory description for the ^5*0 and ^5*1 mesons. A schematic model is 
used for the probability, Pn,qq, of the valence component of the mesons. 

The contribution of the nonvalence component of the photon wave function appears in 
the time-like region, while the nonvalence component of the pion appears in the space-like 
region. The nonvalence contributions of the photon and pion wave functions, relevant for the 
process under consideration, involve emission/absorption amplitudes, that in principle can 
be calculated from the valence components of the corresponding particles, and a suitable 
kernel. However, since our knowledge of this kernel is poor, we use a constant vertex 
approximation for the emission/absorption amplitudes 
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To simplify our calculation, we take advantage of the smallness of the pion mass, which 
is put to zero. Then, only the "Z-diagram" survives in the space-like region. 

We point out that, for m.,^ = 0, only the instantaneous terms contribute to the pion form 
factor. Therefore, in order to evaluate the pion form factor we need the instantaneous vertex 
functions, which we approximate by the full vertex functions times a constant. 

Only a few parameters define our light-front model: the oscillator strength, the con- 
stituent quark mass, the VM meson masses and widths. We use the experimental width and 
mass for the vector mesons, when available. For the radial excitations above 2.150 GeV we 
use the masses of the theoretical spectrum and a single width as a fitting parameter. It is 
worth noting that the results are not markedly sensitive upon different pion wave functions, 
like the asymptotic wave function and the full- model one of Ref. jsil. This could be as- 
cribed to the strong pion binding, that makes the pion wave function similar to its PQCD 
asymptotic limit j55|. 

In the space-like region, the pion electromagnetic form factor is very well described on 
the whole experimentally-explored interval, i.e. up to = — 10 (GeV/c)'^. In the time-like 
region, we find a general agreement up to 10 {GeV/cY, except near the experimental dip at 

2 {GeV/c)\ 

Our model can be straightforwardly improved in many respects. For instance: i) more 
realistic VM wave functions can be used, as the ones of Ref. Q, 

that take into account, 

e.g., the D-state nature of some of the VM resonances, as p(1700); ii) the introduction of 
both a dynamical mixing of isospin states and the contribution of meson. 

Other improvements, like taking care of the non vanishing pion mass, or considering a 
more realistic model for the instantaneous vertexes and for the emission/absorption of a 
pion by a quark, are highly non trivial. 

In summary, our work appears an encouraging step forward in achieving a detailed in- 
vestigation of important issues, as the light-quark content of the photon valence light-front 
wave function, through the analysis of the pion electromagnetic form factor in the time-like 
region. The peculiar feature represented by the smallness of the pion mass is the key point 
to accomplish such an investigation. 
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APPENDIX A: VECTOR MESON DECAY CONSTANT 



The vector meson decay constant, fvn, of the nth state of the vector meson is defined as 



E^v^/y,„=(O|g(O)7^g(O)|0„,A) 



(Al) 



where are the VM polarization vectors and \(f)n,\) is the VM state. 

Let us begin with the four-dimensional representation of the decay amplitude in terms 
of the Bethe-Salpeter vertex of the vector meson, and use the plus component of Eq. ()A1|) 
and A = 2; in the rest frame of the vector meson, where = [M„, 0] and = 1: 



fvn 
Tr 



Nr 



— I 



dk dk^dk. 1 



4 (27r)4 J fc+ {P+ - k+) {k~ - k-, + ^) (P, 
i^-pn + m) 7+ + m) L(P„) ■ [Vn{k, k - P„)], 



k- 



X 



■k+- 



(A2) 



A factor of \/2 enters in the denominator of Eq. ()A2j) from the normalization of the neutral 



vector meson, i.e., (uu-dd)) 1^/2. Then integrating over k , with the assumptions on the 



VM vertex function already presented at the beginning of Sect. Ill, and taking advantage of 
the identification of Eqs. ()76p and ()77|). one arrives at a three-dimensional formula for the 
decay constant where the valence component of the vector meson wave function appears: 

A^, /-A/^ dk+ c/k I 



fvn 

Tr 



4 (27r)3 Jq k+ (M„ - k+) 



^„(A;+,k^;M„,Ox) x 

i^kon (Pn k^on) 



Mo{k+,kr,P+,Pn±) + 2m^ 
Evaluating the trace in Eq. ()A3|) the final expression of the decay constant is: 

Nr r^^" dk+ dk I 



(A3) 



fvn ■ 



k^;M„,Oj 



87r3 Jo k+ (M„ - k+) 

kon ' i^Pn k)ofi (Pji k^Q^ z k -\- kgYi^z ^ ) 

i^kon (Pfi k)on)2 



+m (2fc+ - M, 



Mn + 2m 



(A4) 



where = {k]_ + m'^)/{x{l - x)) with x = fc+/M„. 

From the vector meson decay constant one gets the decay width to e~^e~ as 



Vn 



M3 



(A5) 



where a is the fine structure constant. 
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APPENDIX B: CURRENT CONSERVATION 



Let us define the four quantities 



M2 M2 



where, see Eq. |TBj) of Sec. II, 

[\>(A;,A;-g) 
One can immediately verify that 



7' 



Mo + 2m 



on 



q-Vn = 



Since the vector meson propagator is given by 



-9' + ^ 
M2 



a possible conserved photon- (gg) dressed vertex can be defined as follows : 

n 

with 

fvn 



-9^ + 



<fqu_ 

M2 



V2 VtiAn{k,q) 



KAnik^q) 

fvn 



q'-Ml + iMnTniq' 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



Indeed it is straightforward to show that q ■ J{k, q) = 0, since one has g ■ V„ = 0. 

For each term J'^, let us consider the reference frame where q^ = Mn > and q± = 
( see Sect. VII for the possibility to use different reference frame for different terms in the 
sum). In this reference frame one has q~ = q^ /Mn and then 



q-Vn = - [q-vl + q^v;^ 



VI q^/Mn + MnVn 



(B7) 



Therefore one obtains 



q-Pn 



My I " 



q^Ml 



(B8) 



(B9) 
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and in conclusion we have 



V2 Vn,z An{k,q) 



Vn 



If in Eq. fIBljl the quantity V'^ is replaced by 

v;^(fc, k - P„) : 



r- 



(BIO) 



(Bll) 



Mn + 2m 

as defined in Eq. (jl4p . then the current J^{k, q) is a four vector. 

Let us note that, if in Eq. (jBlOj) the momentum component of the VM vertex function 
A„(fc, q) is taken at the VM pole, then the plus component of the current J^{k, q) coincides 
with the one used in our calculations (see Eq. (jll7|) ). 



APPENDIX C: SUBTRACTION OF THE BARE TERM IN THE PHOTON LF 
VERTEX 



Let us calculate explicitly the contribution to the current of a 7^ bare term in the case of 
massless pions, in coUinear kinematics (q_L = 0, Pttj, = 0, Ps-_l = 0). As discussed at length 
in Section VII and VIII, only the instantaneous terms can give a contribution in the limit 
of a vanishing pion mass. Therefore one has : 



5 j+ (X 1 d^k A^{k - P^, P^) K{k, P^) Tr 



5 7 



+ m\ 7 [f + mj 



X 



(CI) 



Performing the trace, the integration on k~ in Eq. ()C1|) and by using the identifications in 
Eqs. dni), (|Z3) and (jHl), o one has: 



k+ {k+ - q+) 



(C2) 



5 j+ oc dk+ J dk^ ^^^^^ _ _ j^^. ^ 

^^(A;+, k^; P+, 0^) [ml - M^{k^, kr, P^, 0^)]/Pf , 

where P^ = g"*". Using the fraction x = k'^/q'^ in Eq. ()C2|) . one obtains : 

^ f'dx r Mk\ k^; P+, 0^) K - k^; P^, o±)] 

^~ Jo~ J"^^^ [g2-M2(A;+,k^;g+,0^)+ze] ^^^^ 
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The bare term of the current violates the current conservation. Indeed, the matrix element 
of the minus component of the bare term is 

dj oc dk J dk^ ^^^^^ _ _ ^^^^^^^ ^ X 

k^; P+, Ox) [ml - M^{k+, kr, P/, Ox)]/P/ 



lo X J ^ " g+ [g2-M2(A;+,kx;g+,0x) +26] 
Mk^, kx; P^, 0±) [ml - M',{k+, kr, P/, Ox)]/P+ . (C4) 



Then one has 



S q + S j oc 



_ 2 f'dx r P. V;^(fc+,kx;P+,Ox) K - Mg(fc+, k^; P+, 0^)] 
^ io X i ^ [g2-M2(A;+,kx;g+,0x) +26] P/ 

'dx f,,^ ,,2n.+ u .„+ n ^ ^.^.(A;+,kx;P+,Ox) K - M2(A;+, k^; P+, Oj 



+ / ^ /(ikxMo2(A;+,kx;g+,0 
Jo X J 



lo X J " , , [g2-M2(fc+,kx;g+,0x)+2e] P/ 

= _ f'^ Ui,^ P.^^(A;+,kx;P+,Ox) K-Mo^(fc+,kx;P+,Ox)] _^ ^ ^ ^^^^ 

Furthermore, in the present model the pion wave function at large momentum has the 
asymptotic form given by Eq. p45p . which decays as l/|k_Lp. Therefore the matrix element 
of the plus component of the bare term of the electromagnetic current is ultraviolet divergent. 
Therefore in the model of the present paper, which deals with a massless pion, we give out 
the bare term, which gives an unphysical result. 



APPENDIX D: NORMALIZATION OF THE VECTOR MESON WAVE FUNC- 
TION 

The light-front wave function of vector mesons includes the relativistic spin part, as we 
have written in Eq. (jZH)), so one has to consider the whole structure of the wave function 
to obtain the normalization of the valence component of the state. We normalize the wave 
function using the good component of the vector current, imposing that 

a|J + (O)|0„, a) =2MnPg-g,n, (Dl) 

where it appears the probability of the valence component in the vector meson, Pqq-n, which 
is estimated in the next Appendix E. 
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The matrix element of the good component of the current in impulse approximation 
is represented by a Feynman triangle diagram. After the integration over the light-front 
energy, performed disregarding the singularities of the VM vertex in the k~ complex-plane, 
one obtains the contribution of the valence wave-function to the normalization: 

1 iV, l-Mn dk + 

2Wn 16^ io (fc+)2(M„ - k+) 
In Eq. f)D2|) the quantity JV^k^, k_|_) is the following trace 



1 r^^" (]k^ r 



Af{k+,kj 



Tr 



U - Pn)on + m] [/a + {kon ' {Pn - k)on) ' tx Hs{Mo)] (Ikon + m)7+(|^o„ + m) X 



[/a + (kon - (Pn - k)on) " CA ^^(Mo)]] (D3) 



where Hs{Mq) is defined in Eq. 

To evaluate the normalization, we choose the polarization in the transverse direction. 



which is free of the pair term contribution in the limit of zero momentum transfer 
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APPENDIX E: ESTIMATE OF THE PROBABILITY OF THE VALENCE COM- 
PONENT 

In this Appendix we construct in the Fock space of constituent qq pairs a schematic 
model for the light-front square mass operator, M^, which allows one to roughly estimate 
the probability of the valence component in the nth excited vector meson state, Pqq-n- We 
are looking for a square mass operator in the Fock space with a spectrum where the mass 
of the nthth vector meson grows with y/n, as occurs experimentally HQ- 
Let us denote by z > the number of qq pairs and by the noninteracting Fock state 
with i pairs. 

Let us suppose that the free mass operator, Mq, is additive in the number of pairs and 
therefore that the noninteracting Fock-state \i)o is eigenfunction of the free squared mass 
operator with eigenvalue z^: 

M^\i)o = a^ \i)o . (El) 

where a is the energy of a free qq state. 

We suppose that the interaction, Mj, in the squared mass operator (M^ = Mq + Mf): 
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i) has constant diagonal matrix elements 

o(^|M|K)o = ^-c, (E2) 

ii) mixes the state \i)o with the states \{i — l))o and \{i + l))o, and is attractive and 
constant: 

o((2 + l)\Mf\z)o = o{i^ - 1)|M/ K)o = , (E3) 

while the other matrix elements of Mj are supposed to be zero. 
The eigenvalue equation for the squared mass operator is 

Mo' 1^) + \n) = Ml \n) , (E4) 

where the nth excited state of the meson has mass M„. The VM wave function for the nth 
excited state in the Fock space is given by 

l'^) = H K)o ) (E5) 

i > 1 

and the amplitudes a„^i are normalized as follows : 

E l«n,ir=l . (E6) 
i > 1 

In the above sums one has i > 1, since the vector mesons have quantum numbers different 
from the vacuum and then o = 0. 

Introducing the interaction defined by Eqs. ()E2l IE3|) into Eq. (IE4|) and projecting the 
eigenvalue equation in the Fock-space state basis {|i)o}, one has: 



a I an,i = [M^ + c j arr^i . (E7) 



If we define x = A ■ i, with A > 0, then Eq. ()E7|) can be rewritten as follows 

an,i = — [M^ + c) an,i (E8) 

Then defining aS/A"^ = Q/A and going to the continuous limit one gets: 



x'^ttnix) = A„ a„(x) , (E9) 
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with the boundary condition a„(0) = 0, in order to reflect the constraint a„ o = 0. Trivially, 
because of this constraint the eigenvalues and the eigenstates of correspond to the odd 
eigenvalues and eigenstates of the unidimensional harmonic oscillator, namely 



Xr 



n 



2" n\ 



n 

Att. 



1/4 



-e/2 



(ElO) 
(Ell) 



where i?n(0 ^^^^ the Hermite polynomials and 



-X 



The function an(x) in our case has to be normalized as follows 



|a„(x)pc/x 



(E12) 



in order to have the proper correspondence with Eq. ()E6|) . This normalization explains the 
presence of the factor \/2 in the definition of a„(a;) (Eq. (jEllj) ). 

Then, defining a/6 = u/A, from the eigenvalues A„ of Eq. ()E9jl one can obtain the masses 
of the vector mesons, viz 



Mi 



C = UJ (nex + -) - c , 



(E13) 



UJ 1. 

— m + - 

2 ^ 2' 

where only the odd values of n, namely n = 2?2ex + 1, are allowed. The number nex is the 
excitation number 0, 1, 2, 3 ... of the vector mesons, with zero for the meson ground state. 
Then the square mass of the vector mesons is given by: 



Mi 



■ nex + M^ 



(E14) 

where the meson ground state has mass Mg^ = Su/i — c (note that, given u, the constant 
c is fixed). 
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33|, which is the 



At this point we exactly retrieve the experimental spectrum law 
motivation of this simple model. 

The final step is the estimate of the probability of the valence state, namely Pqq,n- From 
a comparison of the discrete and the continuum case of our model, one can associate to the 



probability |a„^j| of the wave function component with i qq pairs, Eq. ()E5|) the quantity 



dx |a„(x)P = 

(i-i)A ' ^ ^' 2"n! 



1/2 .iA 

/ dx\Hr.{0?e 
J{i-1)A 



2"' n\y/Ti J{i-l)S^/2 



4tt 

di \Hrm?e-^ 



(E15) 
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Then the probabihty Pqq,n^^ of the valence component is given by the quantity 

Pg5,n..= |an,i|'~ / dx |a„(x)p= / rf^ l^nlOT e'^ . (E16) 

JO / nlvTT Jo 



Imposing that Pqq-fi, i.e. the valence component probability in the ground state of the 
vector mesons, is about the same as the one found in constituent quark models of the pion 
j2ll |. i.e. it is equal to 0.77, we obtain uj^6 = 2.94. 

Alternatively, let us evaluate the average number of qq pairs in the vector meson: 



V< x2 > 2/1 2 / 3 
Then we can estimate the probability for the lowest Fock component to be roughly given by 

Pqq;n,^ = / • (E18) 

2 \/2nex + 



2 



If the value cja^ = 2 is used in this last estimate, the probabilities considered in Ref. 
are obtained. With this choice, the valence component probability in the ground state of 
the vector meson is equal to ~ 0.8 
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TABLE I: Known vector meson masses, M„, and widths, r„, used in the model. The corresponding 
decay widths into e^e^ pairs, calculated with the VM valence probabilities Pqqyn obtained in 
Appendix E (see also Table II) and the oscillator strength uj = 1.556 GeV^, are compared with the 
experimental values from J34]. ( See text for details) 



Meson 


Mn (MeV) 


M-P (MeV) [34] 


r„ (MeV) 


r-P (MeV) M 




-e- (KeV) 


r;;:^. (KeV) [34] 


p(770) 


770 


775.8 ± 0.5 


146.4 


146.4 ± 1.5 




6.98 


7.02 ± 0.11 


p(1450) 


1497 [42] 


1465.0 ± 25.0 


226 [42] 


400 ± 60 




1.04 


1.47 ± 0.4 


p(1700) 


1720 


1720.0 ± 20.0 


220 


250 ± 100 




0.98 


> 0.23 ± 0.1 


p(2150) 


2149 


2149.0 ± 17 


230 \50] 


363 ± 50 




0.65 
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TABLE II: The vector-meson valence probabilities Pqq-n ^or the first 10 resonances. 



n 





1 


2 


3 


4 


5 


6 


7 


8 


9 


Pqq;n 


0.77 


0.31 


0.29 


0.27 


0.22 


0.18 


0.18 


0.18 


0.17 


0.16 
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TABLE III: Vector meson masses, M„, and widths, F^, used in Ref. 2^]. The corresponding 



decay width into e e"*" pairs calculated with the VM valence probabilities Pnq- 



qq;n 



and the 



oscillator strength uj = 1.39 GeV'^, are reported in the sixth column. (^^The value of 180 MeV for 
the width of /?(2150) is the lower bound of the value obtained by Anisovitch et al quoted in [s^.) 



Meson 


Mn (MeV) 


M^^P (MeV) [34] 


r„ (MeV) 


r^^P (MeV) [34] 




e (KeV) 


(KeV) [34] 


p(770) 


750 


775.8 ± 0.5 


149 


146.4 ± 1.5 




6.37 


7.02 ± 0.11 


p(1450) 


1465 


1465.0 ± 25.0 


310 


400 ± 60 




1.61 


1.47 ± 0.40 


p(1700) 


1723 


1720.0 ± 20.0 


240 


250 ± 100 




1.23 


> 0.23 ± 0.1 


p(2150) 


2150 


2149.0 ± 17 


isot 


363 ± 50 




0.78 
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FIG. 1: Covariant amplitude for tt 7* — tt', or 7* — tttt'. The final meson tt' is a pion in the elastic 
case or an antipion in the production process. 
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FIG. 2: Dressed photon-decay amplitude (7* vrvf) with two possible x"*" time (light-front time) 
orderings, represented by diagrams (a) and (b). The two dashed vertical lines represent different 
light-front times (the light-front time flows from the right to the left). Diagrams (a) and diagram 
(6) contain the processes q ^ qit and q qir, respectively, represented by a full square. The 
crosses indicate the quark lines which are on shell, after the k~ integration. The dashed circle 
represents the dressed photon vertex(see Fig. 4 for details). 
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FIG. 3: Diagrammatic representation of the space-like elastic form factor of the pion for g"*" > 0. 
The light-front time ordering allows one to single out two-quark and four-quark configurations 
at different light-front times, as indicated by the dashed vertical lines. Diagram (a), where < 
—k'^< P+, represents the contribution of the valence component in the wave function of the initial 
pion. Diagram (6), where < k'^ < q'^, represents the non- valence contribution to the pion form 
factor (pair production process) . Both processes contain the contribution from the dressed photon 
vertex. The full square is the vertex function which describes a pion absorption by a quark. The 
dashed circle represents the dressed photon vertex(see Fig. 4 for details). 
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FIG. 4: Dressed photon vertex. The double-wiggly hnes represent the Green function describing 
the propagation of the vector meson V^. The loop on the right represents the VM decay constant, 
fvn (see Appendix A). 
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FIG. 5: (a) Valence component, \qq), of the wave function of an incoming system (photon or 
meson), (b) Non valence component, \qqqq), of the wave function of the same incoming system. 
The extra qq pair is radiatively emitted by a quark in the valence component 
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FIG. 6: Virtual decay amplitude for a pion emission from a quark {q qn) produced by the 
operator /C, see Eq. (jHU- The analogous diagram for pion absorption by a quark can be easily 
obtained by replacing the final pion leg with an initial pion leg. 
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FIG. 7: Space-like em form factor of the pion, tt 7* — >^ tt', for g+ > 0. Two possible light-front 
time orderings are shown. The first one allows to single out the following processes: (a), where a 
point-like quark-photon interaction occurs, and (b), where the absorption of a qq pair by a quark 
proceeds through the kernel 7i. Diagram (c), where the process 7* qq appears (pair production 
process) with the subsequent absorption of the initial pion by a quark, corresponds to the second 
time ordering. 
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FIG. 8: Instantaneous contributions to the time-like em form factor of a massless pion. The 
instantaneous quark leg is attached to the pion vertex in (a) and to VM vertex in (b). The dashed 
circle represents the dressed photon vertex (see Fig. 4 for details). 
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FIG. 9: Pion electromagnetic form factor as a function of the momentum squared . Results for 
the asymptotic and the full pion wave functions, obtained with wvm = —0.7 (see Sect. IX) and 
the quantities shown in Table I, are indicated by dashed and solid curves, respectively. The thin 
solid line represents the result with wvm = — 1/3 and the parameters of Table III. Experimental 
data are from Ref. [3] (full dots) and Ref. |3] (open squares). 
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FIG. 10: Space-like pion electromagnetic form factor divided by the monopole M(g^) = 1/(1 — 
vs the momentum squared . The soHd curve corresponds to wyM = — 1-5 and the dashed 
hne to WyM = —0.7, all the other quantities are according to Table 1. The experimental data are 
as in Fig. |UJ 
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FIG. 11: Pion electromagnetic form factor as a function of the momentum squared q^. The sohd 
curve corresponds to wvm = —1-5 and the dashed line to wyM = —0.7, all the other quantities 
are according to Table I. The experimental data are as in Fig. |^ 
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